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ROUTINE DERIVATION OF POLYNOMIALS TO 
REPRESENT FUNCTIONS 

R obert W. Schrage 

Esso Standard Oil Company 
Linden, New Jersey 

A problem frequently met in digital computer programming is the suitable rep- 
resentation of given functional relationships. Sometimes an analytic expression for 
the function is available and no particular programming difficulty exists. Often, how- 
ever, the function is defined only as a table of corresponding values or as a curve 
plotted on a chart. When this is the case, there are two methods by which the func- 
tion can be converted into a form usable on a digital computer. Either a table look-up 
procedure can be used or an empirical equation may be derived to represent the func- 
tion. Combinations of both methods can also be used successfully. Most of the single 
variable type of functions usually met in technical work can be adequately represented 
by polynomials of a degree no higher than the fourth. Since we have had occasion to 
use a considerable number of such polynomials in our work, it was desirable to devel- 
op a routine method for deriving them. This method has been programmed on the IBM 
Card- Programmed Electronic Calculator so that only a few minutes are needed to con- 
vert a given chart or table into an equation. 

The mathematical technique of using orthogonal polynomials is the foundation 
of the method. Fortunately, it is not necessary to know a great deal of theory to make 
this application.* For our purposes we have so standardized the calculations that a 
person totally unfamiliar with either orthogonal polynomials or computers can easily 
find a polynomial to fit a function. The standardization comes about chiefly through 
an arbitrary decision to select values of the function at a fixed number of equally 
spaced argument values, and to use these function values to derive the desired poly- 
nomial. 

Let us consider the function 

y = f(x). (l) 

It is possible to select several function values at equally spaced values of x and de- 
rive a polynomial curve which will approximately fit them. This curve may be repre- 
sented by 

y = A Q + A lBl + A 2 g 2 + A 3 g 3 + A 4 g 4 (2) 

* See Jack Sherman's paper "The Use of Orthogonal Polynomials in Curve Fitting 
and Regression Analysis" in the IBM INDUSTRIAL COMPUTATION SEMINAR PRO- 
CEEDINGS, September 1950, for a brief, readable introduction to orthogonal poly- 
nomials. 



where the A's are constants and the g's are first, second, etc. , degree polynomials 



in x. 



If these polynomials are of a special type known as orthogonal polynomials, 
then equation (2) will give a least squares fit of the original function values. This 
fit will correspond to the highest degree polynomial included in equation (2). The 
values of the constants, or A's, are independent of the number of terms in equation 
(2). 

Let us assume that we will deal specifically with seven data points. Then the 
values of the A's up to A . may easily be obtained from the following matrix-by-vector 
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Having the A's, the original function is now represented from equation (2) by any de- 
gree of polynomial up to and including the fourth. 

However it is not convenient to express the g's in terms of x, the original 
independent variable. Rather than transform x to these several g values every time 
we wish to use the equation, it would be preferable to transform equation (2) once and 
for all into a polynomial directly in x. 

We can do this in two steps by further matrix-by -vector multiplications. First, 



we evaluate five constants a. 



a, by the operation 
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These a values are the coefficients of an equation 
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y = a Q + ajZ + a g z +a,z + a.z' 
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In equation (6) x 4 is simply the argument value of the fourth or central point in the 
original seven and m is the constant interval between points. It is frequently con- 
venient to use equation (5) directly because it is easy to evaluate z from equation (6). 
Furthermore, z has the valuable property that it always lies between ± 3 if equation 
(5) is used only in the range of the original seven points. This fact is particularly 
helpful when the polynomial is to be used on a fixed decimal calculator. 



If a floating decimal calculator is to be used, however, then it is useful to 
proceed a step further and calculate the five coefficients b Q , 



b, from 
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The b values calculated in equation (7) are the coefficients of a polynomial 

directly in x: 

2 3 4 

y = b Q + b 1 x + b 2 x + bgX + b 4 x (8) 

There may be times when equation (8) can be used successfully on a fixed decimal 
calculator. This depends, of course, on the scale of x. 

It is important to note that although the coefficients A Q , . . . , A. are indepen- 
dent of the number of terms in equation (2), this is not true for a„, . . . , a 4 or 
b Q , . . . , b 4 . These latter coefficients will be different for different degree polynom- 
ials. However, values for degrees less than the fourth can easily be calculated by 
arbitrarily fixing A., or A. and A„, etc. , as zero. 

The matrices in equations (3) and (4) are constant so long as only seven data 
points are used. If it were desired to use fewer or more points than seven, then new 
matrices would have to be derived. However, similar techniques would be used. The 
matrix used in equation (7) does not really depend upon the number of data points. 



It is merely the means for transforming a polynomial in the form of equation (5) to 
that of equation (8). Therefore, provided x. and m are generalized to mean the type 
of constants in equation (6), the matrix in equation (7) does not vary. 

Since only small matrix-by-vector multiplications are involved in these calcu- 
lations, they may be done conveniently on any medium scale calculator. We have de- 
veloped a program for the CPC using a general purpose floating decimal control panel. 
This program derives coefficients for all the types and degrees of polynomials dis- 
cussed in this paper. It also calculates check values so that the fit of the derived 
polynomials to the original function can be examined. Less than five minutes of cal- 
culator time is required for execution of the program. 

An outline of the procedure followed in the program is given in the Appendix. 
A numerical example is also given, showing the various coefficients derived in a 
typical curve fitting problem. Figure 1 shows graphically the original seven data 
points used in the problem and a curve representing the fitted fourth degree polynom- 
ial. Within the accuracy limitations of the drawing the curve coincides with the ori- 
ginal functional relationship. 

It may be worthwhile to emphasize the value of calculating check points inter- 
mediate between the original data points. Sometimes polynomials may be derived 
which fit the given data points well enough, but which do not adequately represent the 
function between them. This difficulty usually arises when the selected points do not 
adequately describe the curvature of the function. The remedy is either to use a 
procedure based on more points (and possibly higher degree polynomials) or to divide 
the function into two or more intervals to be represented by different polynomials. 
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APPENDIX 



Program Outline 



Reading of problem data. Seven 
function values, taken at equally 
spaced values of the argument, are 
read. The argument value corres- 
ponding to the first function value 
and the increment in argument values 
are also needed. 



Calculation of the orthogonal poly- 
nomial coefficients used in equation 
(2). This is done by the matrix-by- 
vector multiplication given in 
equation (3). 

Calculation of coefficients for the 
type of polynomial defined in 
equations (5) and (6). Complete 
sets of coefficients for all degrees 
of polynomials up to the fourth are 
calculated. This is done by the 
matrix-by-vector multiplication 
given in equation (4). 



Numerical Example 
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3rd degree: 
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4th degree: 

If, 3 

a, = -2 

a; = 1 

ar = -1 



,009784 
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4. Calculation of coefficients for the 
final form of polynomial as given 
in equation (8). The coefficients 
are obtained by the matrix-by-vector 
multiplication given in equation (7). 



1st degree: 



h0 = 
b l = 



2nd degree: 
b = -1. 
b" = 7. 

b 2= " 7 - 
3rd degree: 
b = -1. 
bV = 1. 
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b 3 = 9 ' 
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Program Outline 



Numerical Example 



4th degree: 
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Calculation of check values. These 
are obtained by substitution of the 
appropriate argument values in each 
of the several degrees of poly- 
nomials derived in step 4. * Since 
these last polynomials depend upon 
the correctness of all preceding 
calculations, this is the best place 
to make the check. 

The check values would be the same 
for any form of polynomial - 
equations (2), (5), or (8) - of 
equivalent degree. 



Check values for 4th degree: 
6.194140 xl0-j 
7.317100 xlO 
8.223130 x 10 
8.934540 xlO 
9.471850x10 
9.853850 x 10 
1.009778 x 10 
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1.023136 x 10 
1.014670x10 
9.975520 x 10 
9.726300x10 
9.405970 x 10 
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* To save space, check values for only the fourth degree equations are shown in 
the.numerical example. The intermediate points are at argument values mid- 
way between original data points. 
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AN INFORMATION SEARCHING SYSTEM BASED ON 
THE IBM COLLATOR 

E. R. Lancaster 

Department of Mathematics 
Newark College of Engineering 



From each of a group of documents is selected a set of words which more or 
less identifies the contents of the document. We prepare a file of these "key" words 
in such a way that the file can later be searched for an arbitrary set of words and a 
list of documents obtained, each of which contains all of the words in the set. 

The manner in which the file is prepared will be discussed after its operation 
is described. 

Operation of the System 

The file consists of IBM punched cards. Following is a column breakdown of 
the information punched in a typical card prepared for a magazine article. The en- 
tries for the first seven columns will be different for other types of documents, such 
as engineering reports, patent digests, etc. 

Columns Entry 

1-3 Source code. 3-letter code for the name of 

the periodical in which the article appears. 

4-5 Last two digits of the year in which the 

article appeared. 

6-7 Number of the month or week in which the 

article appeared, e.g., 11 for Nov. , 31 for 
July 21, 1954. 

8-21 Key word. There is a card for each key word. 

The placement of words of more than 14 letters 
in these columns is explained later. 

22 Number of card for a given key word. Occasionally 

a title cannot be fitted into columns 23-80, in which 
case a second card is necessary. If only one card 
is required, this column is left blank. Further 
discussion will be given later. 

23-80 Title of article. 

The file of cards is arranged in alphabetical order by key word. All cards 
having the same key word are arranged in numerical order on columns 1-7, 23-24. 
Cards which are identical in columns 1-21, 23-24 are further arranged in numerical 
order on column 22. 
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Suppose we wish to obtain a bibliography on friction clutch materials. First, 
we refer to a list which gives the number of times each key word appears in the file. 
If all three of our key words (friction, clutch, material) appear there, we manually 
remove from the file all cards for these key words. Suppose that of the three decks 
removed the "friction" deck contains the fewest cards. This deck is placed in the 
secondary feed hopper of a collator and either of the other decks in the primary feed 
hopper. Assume we place the "clutch" deck in the primary hopper. The control 
panel is wired in such a way that the primary deck falls in pocket 2, cards in the sec- 
ondary deck which match in columns 1-7, 23-24 with cards in the primary deck fall 
in pocket 3, and the remainder of the secondary deck falls in pocket 4. When this 
process is completed, we will have in pocket 3 the cards for those articles from which 
"friction" and "clutch" were originally selected as key words. 

The "clutch" deck is taken from pocket 2 and returned to the file. The cards 
in pocket 4 are placed temporarily in a special rack. The cards in pocket 3 are 
placed in the secondary hopper, and the "material" deck is placed in the primary 
hopper. The cards in the "friction" deck which match with cards in the "material" 
deck again fall in pocket 3, and these are the cards for those articles from which 
"friction, " "clutch, " and "material" were originally selected as key words. The 
"material" deck is taken from pocket 2 and returned to the file. The cards in pocket 
4 from the "friction" deck are placed temporarily in the special rack. The informa- 
tion in columns 1-7 and 22-80 of the "friction" deck cards in pocket 3 is listed by 
means of an IBM accounting machine. If an accounting machine is not available, a 
list can be made with a typewriter. 

This list is the completed bibliography of articles on friction clutch materials. 
It gives the title of the article, a 3-letter code for the periodical in which it appeared, 
and the date of publication. Although author's names do not appear in the printed bib- 
liography, they can be used as key words in any search. Author's names appear as 
key words in the file for every article covered by the file. In the list, the appearance 
of a 2 just to the left of a title will indicate it is to be read as a continuation of the 
previously listed title. 

Upon completion of the listing, the cards from the "friction" deck which were 
placed in the special rack are merged, by means of the collator, with the cards used 
in the listing. The resulting deck (the "friction" deck) is returned to the file. 

All the cards for this system may be placed in one file. However, when the 
number of cards becomes extremely large, it may be more efficient to have several 
files. This division might be made by classifying the documents covered by the sys- 
tem. A simple example is the division of the cards into mechanical and electrical 
files. Any such division by subject heading would require the placement of cards for 
some documents in more than one file. A decision as to which files to search in the 
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preparation of a particular bibliography would be made after examining the lists to 
determine the number of times each key word appears in the various files. 

Preparation of the System 

In the following discussion we will continue to consider a magazine article as 
a typical document. 

A person, whom we call the analyzer, selects the key words for a given arti- 
cle. He then writes the results of his analysis on special paper which is ruled into 
squares in such a way that the squares are easily counted. 

The analyzer first rewrites the title of the article with as few letters as pos- 
sible, using one square of the paper for each letter and space between words. If the 
revised title has a length greater than 58 squares, letters after the 58th square are 
placed on the second line. 

After the title the analyzer writes on another line the information which will 
appear in columns 1-7 of the card. He then scans the article for key words, which 
he underlines and writes on successive lines of the special paper after changing to a 
basic form those not already in basic form. For example, "experimental" would be 
changed to "experiment," "verification" to "verify, " "properties" to "property." 
As the work of preparing the file progresses, a dictionary of basic forms is compiled. 
This dictionary also includes a list of synonyms and other related words to help in the 
selection of key words for file searches. For example, the synonyms "computer" 
and "calculator" would be placed in the same word group in the dictionary and only 
one would ever appear in the file. 

Key words exceeding 14 letters in length are decreased to 14 by marking out 
letters from the right in accordance with the following scheme until the length is 14: 

Number of letters Mark out every 

8th 
5th 
4th 
3rd 
2nd 

Thus "thermoelectrodynamometer" becomes therolcrdnmmtr. " 

The analyzer places the completed sheet in one of two stacks, depending on 

whether one or two lines were required for the title. When the stacks have become 

sizable, they are given to a key-punch operator. 

The key-punch operator first prepares a deck of cards for the stack of sheets 

with one-line titles. From each sheet are prepared one master card for that sheet and 

a detail card for every key word on the sheet. On the master card the indicated letters 



15 to 16 


17 


18 


19 to 21 


22 to 28 
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and digits are punched in columns 1-7, an X-punch in column 22, and the title in col- 
umns 23-80. The key words are punched, one per detail card, in columns 8-21. The 
resulting deck is then placed in a reproducer and the information in each master card, 
with the exception of the X-punch in column 22, is reproduced in the same columns of 
the detail cards following it. The master cards are removed with a sorter and the re- 
maining cards, called deck A, are stored until work is completed on the sheets with 

two-line titles. 

The key punch operator must go twice through the stack of sheets with two-line 
titles. The first time through, she punches a master card for each sheet and a de- 
tail card for every key word on the sheet. On the master card, columns 1-7 are pun- 
ched, an X-punch is placed in column 22, and the first line of the title is punched in 
columns 23-80. Following each master card, a detail card is punched for each key 
word from the article, as already described for one-line titles. The deck resulting 
after this first run through the sheets is called deck B. Deck B is now run through a 
reproducer, and the punches in columns 1-7 and 23-80 of each master card are repro- 
duced in all the detail cards between it and the next master card. Following this, Deck 
B is again run through the reproducer and a new deck, Deck C, is made which is a 
duplicate of Deck B in columns 1-22 and blank in columns 23-80. After the master 
cards are removed with a sorter, Deck B is temporarily stored. 

The key punch operator now goes through the sheets with two-line titles a sec- 
ond time, punching one master card but no detail cards for each sheet. These cards 
have been prepunched with an X-punch and a 2-punch in column 22. The operator 
punches only the second line of the title on the master card for that title. With a col- 
lator these master cards are inserted in Deck C immediately following the master 
cards already in Deck C. (If the collator is Type 77, the cards should be fed 12's 
edge first.) 

Deck C is now placed in a reproducer, and the punches in columns 23-80 of each 
master card are reproduced in each detail card between it and the next master card. 
The master cards are then removed from the deck with a sorter. 

Decks A, B, C and the file itself are next sorted into the proper order on the 
following columns in the order given: 22, 7-1, 21-8. 

The single resulting deck is the new file. This deck is placed in an IBM alpha- 
numeric accounting machine, and the words in columns 8-21 are listed without repet- 
ition in one column and the number of times the word occurs in another column. If the 
accounting machine is connected to a reproducer, each word and the number of times 
it occurs can be punched in cards, one word per card. This facilitates further listing. 

If the deck obtained from the reproducer has a color different from the main 
deck, it can serve another useful purpose. We can, by means of a collator, insert 
these cards into the main deck in such a way that cards with the same word in columns 
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8-21 are easily identifiable. This will be helpful in the manual extraction of card 
groups from the file. 

A Useful Auxiliary File 

Deck D consists of 792 cards. The first card in this deck has a single hole 
punched in the 12's row of column 15, the second card has a single punch in the ll's 
row of column 15, the 13th card has a single punch in the 12's row of column 16, the 
32nd card has a single punch in the 5's row of column 17, etc. The 792nd card has a 
single punch in the 9's row of column 80. 

An auxiliary deck is prepared by means of a collator and a reproducer. First, 
the main file of cards is placed in numerical order on columns 1-7. This deck is 
divided into 792 parts, each part having approximately the same number of cards. In 
front of the first part is placed the first card of deck D, in front of the second part is 
placed the second card of deck D, etc. The resulting deck is called deck E. 

Deck F is obtained from deck E in such a way that a given card in deck F has 
in columns 1-14 a key word from a card of the main file plus the single punch from the 
card of deck D which precedes the part of deck E containing the card from the main 
file. Deck F is arranged in alphabetical order on columns 1-14 to obtain deck G. 
All the cards for a given word in deck G are reproduced on a single card. The resul- 
ting deck is the auxiliary deck. It has one and only one card for each key word in the 
main deck. 

The use of the auxiliary deck will be explained by means of an example. Sup- 
pose, as in our previous example, we desire a bibliography on friction clutch mater- 
ials. Before searching the main file as previously explained, we extract from the 
auxiliary file the three cards punched respectively with the words "friction, " 
"clutch, " and "material. " The three cards are placed together in such a way that a 
punch in columns 15-80 in the same position in all three cards will be apparent. The 
absence of such a punch is proof that a search of the main file will not reveal an item 
for the bibliography. 
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THE STANDARD SCORE METHOD OF COMPUTING AN INTER-CORRELATION 
MATRIX ON THE 604 ELECTRONIC CALCULATING PUNCH * 

K. Warner Schaie 
University of Washington 

One of the constantly recurring problems in the analysis of social science 
data is the computation of product-moment correlation matrices. If the correlation 
data are to be used for further analysis, such as factor analysis or regression, it 
is necessary to compute the complete matrix, i.e. , the correlation of each variable 
with every other variable. It is clear, then, that when we consider more than four 
or five variables, manual methods of computation become enormously time-consum- 
ing and often prohibitive, even when the number of cases is relatively small. 

If we deal with a thousand or more cases, or a hundred or more variables, 
it will be economical to program and use one of the larger data processing machines 
such as the 701, 704, 705, or perhaps the 650. In the social sciences, however, the 
typical problem is not of such a magnitude as to warrant use of these larger machines. 
Here, in a typical problem we have about 50 to 500 cases, and we are interested per- 
haps in the correlations of 20 to 50 variables. 

Where it is impracticable to use the larger equipment, the use of the 604 
Electronic Calculating Punch holds considerable promise, particularly since more 
than 2, 000 of these machines are available throughout the country, their programming 
is relatively simple, and their cost of operation is fairly reasonable. 

This paper describes a method by which it is possible to compute 40th order 
and larger correlation matrices at an average rate, for a sample of 100 cases, of 
approximately 100 correlation coefficients per hour, depending upon the skill of the 
operator. This method will result in the complete computation of the inter-correla- 
tion table. Each correlation will be punched to three decimals on a separate card, 
each of which will also be punched automatically with the correct row and column de- 
signation. The check for the computation consists of the independent recalculation of 
the infra-diagonal values and of the calculation of the diagonal. 

THE STANDARD SCORE METHOD 
Let us consider the formula for the Pearson Product Moment correlation 

* This research was part of an inter-disciplinary study supported by a grant from the 
Agnes H. Anderson faculty research fund to Dr. Charles R. Strother, Chairman, 
Committee on Gerontology, University of Washington, Seattle. The writer wishes 
to acknowledge his indebtedness to the staff of the IBM Seattle sales and service bur- 
eaus for their assistance during the machine -testing of the method described in this 
study. 
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coefficient to see how it can be written in steps simple enough for efficient machine 
calculation. To begin, 

xy " Na x a y ' {1) 

where x is of the form (X-X) ancl y is of the form (Y-T). Equation (1) indicates that 
the correlation coefficient between two variables is obtained by summing the cross- 
products of the deviations of each pair of raw scores from their respective means and 
dividing the obtained sum by the product of the standard deviations of the two varia- 
bles times a constant. 

It is easily seen that the storage capacity of the calculator (four three-digit 
units, five five-digit units) as well as the limitations on calculating capacity preclude 
the solution of this equation in a single operation. A solution well within machine cap- 
acity becomes feasible, however, by considering an alternate form of equation (1). 



Let z be a standard score of the form 



z__ = 



X-X 



X CT x 



(2) 



and let z = Y Y . (3) 

y (j 

y 

Then, we can insert (2) and (3) in equation (1) as follows: 

S z z 
i x y 

r xv = i (4 > 

^ Na or 
z z 
x y 

But the standard deviation of any distribution of z -scores equals unity and the denom- 
inator term accordingly reduces to the constant N. Furthermore, since the mean of 
a distribution of z-scores equals zero, all values entering the sub-products of the num- 
erator term will be distributed about zero. 

Since the 604 will recognize a negative number (if it is identified by an X-punch 
over the unit position) and will automatically adjust the sign of products of two nega- 
tive numbers or of a negative and a positive number, the above development suggests 
ways of providing sufficient storage space for the simultaneous development of sev- 
eral products over a substantial number of cases, providing that raw scores are first 
reduced to standard score form. 

The number of cases which can be handled by this method is determined by 
the number of digits to which standard scores are to be computed and, also, by the 
number of digits to which the sum of products is to be accumulated. The number of 
correlations which can be computed at one time is limited by the number of storage 
units available on the calculator. If the maximum of four coefficients at one time is 
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to be computed, standard scores must be expressed in not more than three digits. 
The maximum number of cases which can be handled can be ascertained by equation 

l(T a 

where the superscript a is the number of digits to which standard scores are ex- 
pressed, the superscript b is the number of decimals to be dropped from each pro- 
duct of two standard scores, and N is an integer, any decimal remainder being drop- 
ped. For example, if we plan to compute standard scores to one digit and if we plan 
to drop no decimals, then a = 1 and b = 0, and inserting these values in (5) gives 

N= MM™ =1000.00 . 
100 

Dropping decimals, we would know that in this situation we could handle up to 999 
cases. 

By full use of punch and calculator selectors available on the 604, it is fur- 
ther possible to compute four sets of four correlations in one machine run. This oper- 
ation involves the reproduction of cards punched with the standard scores, but makes 
possible a considerable saving in machine time. In outlining this method we shall use 
an example of a 40 x 40 correlation matrix for an N of 50 cases (standard scores com- 
puted to one decimal), which has actually been computed. A simple modification will 
also be described by which our method can be expanded to handle up to an 80 x 80 ma- 
trix, computed from two basic IBM cards per case. Further modifications are pos- 
sible, which permit the handling of a 200 x 200 matrix from five basic cards. The 
maximum utility of our method, however, seems to find an upper limit in an 80th or- 
der matrix. 

For convenience in presentation, the required procedure will be divided into 
three parts: (1) conversion of raw scores to standard scores, (2) reproducing the 
standard score deck and preparing the operating deck, and (3) calculation of correla- 
tion coefficients from the operating deck. A further section describes the necessary 
modifications for calculating an 80 x 80 matrix from two basic cards, while the final 
section will present a rule of thumb for determining the machine time required to com- 
pute a correlation matrix of any given dimension. 

1. Conversion of Raw Scores to Deviation Scores 
Data Preparation 

This method requires that means and standard deviations be known for each 
variable to be correlated. These values may be either manually computed or obtained 
in a separate operation on the 604 or other equipment which will accumulate sums of 
squares. 
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Raw scores are punched on one or more basic NX cards for each subject. 
Ideally, raw scores should have no more than two digits. In this case, all 40 scores 
for a given subject can be punched on one basic card. It is possible to accommodate 
a few three-digit raw scores (provided the highest digit is identical for all scores) 
by using an X overpunch. X must never be punched in the unit position, however, as 
this position is required to identify negative scores. 

X-21 cards are punched, one for each variable to be correlated. Columns 
1-2 are used to identify the variable; columns 71-75 are assigned to the mean; and 
columns 76-80 to the standard deviation. 

A set of blank cards is punched with a 12-punch in column 61 to receive the 
deviation scores. One 12-61 card is required for each subject plus two 12-61 cards 
required for programming operations. 

12-61 cards are sorted in with the NX cards so that one 12-61 card follows 
each NX card. The following card sequence is required for the conversion operation: 

12-61 (blank) 

X21 (different for each run) 

12-61 (blank) 

"KTV 

(one pair per subject) 
12-61 

Conversion Operation 

Figure 1 shows the planning chart for this computation and indicates that an 
internal check is provided by reverse calculation and zero check before the results 
are permitted to punch. Figure 2 shows the appropriate wiring diagrams. Wiring 
condition (1) is a sample of emitting the third digit from an X overpunch. 

One run is required for each variable to be correlated. Program advance re- 
quires only wiring (2) and (3) to be moved two columns on first reading and punching 
exits, respectively, after each run. 

The same deck of cards is used for each run except that the X21 card is 
changed after each run to feed in the mean and standard deviation for the variable 
which is being converted. 

Note that we suppress punching on NX cards. If we failed to do this, our 
basic deck would be destroyed during the operation. Punching is done out of general 
storage rather than the counter to obviate wiring which would otherwise be required to 
suppress calculation on the NX cards read cycle. 
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2. Preparing the Operating Deck 

In order to minimize program advances, we now reproduce the 12-61 deck 
three times. The 12-punch on these reproduced decks must be transferred, however, 
to 12-65, 12-67, and 12-69, respectively. The reason for this change will become 
apparent later on. Next, we need a set of cards, one for each variable, each having 
the variable designation punched in columns 1-2 of the card, and an X-punch in col- 
umns 61 and 79. For smoothness of operation it might be worthwhile to duplicate the 
complete set so far described, so that one set can be arranged in sequence while the 
other is feeding through the calculating punch. Note that 12, rather than X-punches, 
must be used if negative numbers are to be handled. 

Next, we need a set of blank cards to receive the correlation coefficients 
and their variable identifications. Since we will calculate four correlations at one 
time, we need four different sets to program the separate conditions. We shall de- 
signate these cards XI, X3, X5, and X7, respectively. The number of these cards 
needed is obtained by dividing the total number of correlations to be computed (both 
original and reverse calculation) by four. In our example we have 40 x 40 or 1600 
correlations, which means that we need 400 each of the XI, X3, X5, and X7 cards. 

The card sequence for the final computation will be as follows: 

X7 (blank) 

12-61 (one for each subject) 



X79 

XI 

X3 

X5 

X7 

12-65 

X79 

XI 

X3 

X5 

X7 



(different for each run) 



(one for each subject) 
(different for each run) 



12-67 

X79 

XI 

X3 

X5 

X7 

12-69 

X79 

XI 

X3 

X5 

X7 



(one for each subject) 
(different for each run) 



(one for each subject) 
(different for each run) 



A new blank set of XI, X3, X5, and X7 cards is used for each run. For a 
sample of 50 subjects this would imply an operating deck of 221 cards. Of these, 
only 20 cards would be changed on each run, i.e. , we would need to insert new cards 
to replace the ones on which correlations have been punched and also to insert the 
appropriate variable designations given by the X79 cards. For convenience in making 
these changes, the X79 and X1-X7 cards should be punched on distinctive colors. 
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3. Computing the Correlation Matrix 

Figure 3 gives the basic planning chart for this operation. Standard scores 
are read in through the four three -position storage units and the MQ unit, and the 
sums of products are accumulated in the four five-position storage units. The capa- 
city of these units limits the size of the sample that can be manipulated in one opera- 
tion. 

The XI -X7 cards are used to select the sum of products to be divided by the 
sample size, which will result in the correlation coefficient. N (the sample size) is 
emitted from the digit emitter during division. If this divisor has more than one digit, 
it will be necessary to add the required number of programs to develop the divisor in 
the counter and transfer to a free three-position storage unit. If the divisor has more 
than three digits, it will be necessary to program division in two separate program 
steps. 

Figure 4 gives the wiring diagram for this operation. Note that punching is 
suppressed on all but XI -X7 cards. This design provides for the computation of 16 
correlations in one machine run. This is accomplished by successively computing the 
correlation between the set of four variables A, B, C, and D and the variables E, F, 
G, and H, which are read into the MQ unit successively by wiring through punch se- 
lectors. This is the reason why we reproduce the original standard score set (12-61) 
with the classifications 12-65, 12-67, and 12-69, respectively. These 12-punches 
select the multiplier variable with which the four first variables now in storage are to 
be multiplied. 

Variable designations are provided by reading in the row designation from the 
X-79 cards at Second Reading and punching directly from Second Reading on XI, X3, 
X5, and X7. On the first run we must eliminate reading from columns 1 and 2. On the 
NX card this is done by punching an X-79 into NX01 and wiring FS1, using punch se- 
lector 1 transferred, picked up by X-79. This wiring may be discarded after the first 
run. Column designations are punched from digit emitters which are selected through 
punch selectors picked up by XI -X7 cards to punch correct designation. 

Program Advance Schedule 

In figure 4, solid lines are wiring which will remain throughout the operation; 
dashed lines indicate wires which must be changed for each set-up; and dotted lines 
indicate wires which must be changed for each run. In our example, 16 wires must 
be moved for each set-up, while 8 wires are moved for each run. 

The number of set-ups required is determined by dividing the number of 
variables to be correlated by four. Each set-up again has as many runs as there are 
set-ups. For the 40 x 40 matrix, then, we need 10 set-ups with 10 runs each, or a 
total of 100 machine runs of 421 cards. The program advances involve 9 major and 90 
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minor changes, each of which, however, should not take more than a minute or two. 
Figure 5 shows a sample program-advance schedule for one set-up. 

MODIFICATION FOR COMPUTING AN 80 x 80 INTERCORR ELATION 

MATRIX 

The method described for the 40 x 40 table can be used for the larger matrix 
except that one additional punch selector must be wired on all set-ups where the mul- 
tiplicand and multiplier are read from two different cards. It is also necessary to 
wire the digit selector through a pilot selector to get the correct variable designation. 
In the 80 -variable case, this will involve changing 16 instead of 8 wires for each set- 
up. 

By using the formula described earlier, we find that in the 80-variable case 
we require 20 set-ups of 20 runs each, or a total of 400 runs. Half these runs will 
compute the intercorrelations on one card, and thus there will be 221 card runs for 
a sample of 50 cases. The other half of the operation will involve runs with informa- 
tion read from two basic cards, and thus there will be 421 card runs. 

The number of runs can be halved if we introduce some relatively simple 
wiring modifications. The principle involved in this change is that we will success- 
ively compute the correlation between a set of four multiplicands on card 1 and a set 
of four multipliers, first on card 1 and then on card 2. On the next run we simply re- 
verse this operation and calculate the correlations between a set of four multiplicands 
on card 2 and sets of four multipliers first on card 2 and then on card 1 . The change 
between the two halves of this operation merely involves transposing the wiring from 
the normal and transferred exits of the pilot selectors, determining whether card 1 or 
2 is to be read. 

A special function punch is needed to activate the pilot selectors. It is nec- 
essary, therefore, to emit a 12-61 punch into all cards 1, and a 12-63 punch into all 
cards 2. In setting up the operating deck for the two-card operation, it should be re- 
membered that we need first a deck of the same composition as described in the pre- 
ceding section, and then another deck obtained by collating decks composed of cards 
1 and cards 2, respectively. 

With this modification we will require 10 set-ups of 20 runs each, where 
each run will involve 641 cards and where 32 correlations will be computed in each 
run. Figure 6 shows the additional wiring required to control the 80-variable condi- 
tion. Program-advance schedules can easily be designed by the appropriate expan- 
sion of figure 5. 
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PROGRAM ADVANCE SCHEDULE (Correlation from standard scores) IBM 604 



SET-UP: # 4 


PS 5T to DS 13 


FS 1 to RI 25 &26 


GS 1 to RI 29 & 30 


PS 6T to DS 14 




PS 7T to DS 15 


FS 3 to RI 27 & 28 


GS 3 to RI 31 & 32 


PS 8T to DS 16 


1st run 

PS 2N to col. spl. # 1 & RI 2 
PS2Tto " " #2 &RI4 
PS 3Tto " " #3&RI6 
PS4Tto " " #4&RI8 


6th run 

PS 2N to RI 41 & 42 
PS 2T to RI 43 & 44 
PS 3T to RI 45 & 46 
PS 4T to RI 47 & 48 


2nd run 

PS 2N to RI 9 & 10 
PS2TtoRI 11 &12 
PS 3T to RI 13 & 14 
PS4TtoRI 15&16 


7th run 

PS 2N to RI 49 & 50 
PS 2T to RI 51 & 52 
PS 3T to RI 53 & 54 
PS 4T to RI 55 & 56 


3rd run 

PS 2N to RI 17 & 18 
PS 2Tto RI19 &20 
PS 3T to RI 21 & 22 
PS 4T to RI 23 & 24 


8th run 

PS 2N to RI 57 & 58 

PS 2T to RI 59 & 60 

PS 3T to col spl. # 5 & RI 62 

PS 4T to RI 63 & 64 


4th run 

PS 2Nto RI25&26 
PS2TtoRI27&28 
PS 3T to RI 29 & 30 
PS 4T to RI 31 & 32 


9 th run 

PS 2N to col. spl. # 7 & RI 66 
PS2Tto " " #8&RI68 
PS3Tto " " #9&RI70 
PS 4T to RI 71 & 72 


5 th run 

PS 2N to RI 33 & 34 
PS 2T to RI 35 & 36 
PS 3T to RI 37 & 38 
PS 4T to RI 39 & 40 


10th run 

PS 2N to RI 73 & 74 
PS 2T to RI 75 & 76 
PS 3T to RI 77 & 78 
PS 4T to RI 79 & 80 



Special Wiring on Run # 4 Split-wire 

PS 2N to FS 1 PS 2T to FS 3 PS 3T to GS 1 



PS 4T to GS 3 



Other special wiring conditions: 

Symbols: PS - Punch Selector; FS - Factor Storage Entry; GS - General 



DS - Digit Selector; T - Transferred Side; 
RI - Read in First Reading; N - Normal Side; 
col. spl. - Column Split 



Storage Entry 



Figure 5 
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MACHINE TIME REQUIREMENTS 

The time required for this method is composed of three basic components. 
First, a relatively constant amount of time will be required to wire and test boards; 
the second component is the time required for program advances between runs; and 
the third component is the time used by the machine to run through the required cards. 

In our experience, two hours were required to wire and test the two sets of 
boards. In developing a rule of thumb for estimating calculating time, we shall de- 
fine this time as the constant C. 

The second component may be determined by noting that approximately 1 1/2 
minutes are required for each program advance during the raw-score to standard- 
score operation, and approximately 2 1/2 minutes during the final computation. We 
shall define these constants as A and B, respectively. 

The third component is determined by noting the machine capacity of 100 cards 
per minute. 

From the preceding sections and above information we can write the follow- 
ing formula estimating the required machine time: 

Machine time = n ( 2N 10 + 3 ) + An + r ( 4i [^ 21 ) + Br + C, (6) 

where A, B and C are defined above, N is the number of cases, n is the number of 
variables, and r = n /16, the number of machine runs. 
For our example these values would be: 

MT= 40 (JgJ-) + 60 + 100 (|21) + 250 + 120 

= 41 + 60 + 221 + 250 + 120 
= 692 minutes or 11 1/2 hours 

SUMMARY 

The present paper outlines a method of computing intercorrelation matrices 
on the 604 Electronic Calculating Punch by means of separate conversion of raw scores 
to standard scores. 

Correlation coefficients are punched on individual result cards on which both 
row and column designations are also punched in the same operation. Correlations 
are checked by independent reverse calculation of each coefficient. 

This method appears to be most economical where a large number of varia- 
bles must be correlated for small size samples. Its major utility seems to be restric- 
ted for samples not exceeding 500 and from 20th to 80th order matrices. 
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A METHOD OF CARD PROGRAMMING THE 604 FOR 
GENERAL ENGINEERING CALCULATIONS 
R. Baldur 

Dominion Engineering Works, Ltd. 
Montreal, Canada 



Introduction 



The idea of using the IBM Type 604 Electronic Calculating Punch as a Card- 
Programmed Calculator is not new. However, the methods described previously 
(References 1 and 2) have only a limited application for engineering work, since they 
do not enable the evaluation of the frequently occurring trigonometric functions, and 
they have a limited storage capacity. The following factors influenced the design 
described in this article: 

1. It was required to be able to evaluate the maximum number of the most 
common mathematical functions. 

2. A reasonable storage capacity was necessary to enable easy manual col- 
lation and sorting of cards after each pass. 

3. In order to limit the duplication of the program deck, a method of storage 
by gang- punching had to be avoided. 

4. Double punching had to be reduced to a minimum to simplify the key-punch- 
ing operation. 

Owing to the limitations of the machine, some compromise had to be accepted. 
In particular, the design considered here is only capable of dealing with the four ba- 
sic arithmetic operations, evaluating the square root, and natural and inverse circu- 
lar functions. It is not possible to calculate the hyperbolic, exponential and logarith- 
mic functions, but a table look-up operation is available, which can to some extent 
compensate for that shortcoming. Three storage locations are always available for 
retaining factors or intermediate results. But as many as five can be used if only ad- 
dition, subtraction, multiplication and division is performed, and four when the square 
root is calculated. All punching is normally suppressed, except for the results, when 
four distinct numbers are available. Finally, the design of the cards is such that 
there is only one card which requires double punching in any table look-up operation. 

In order to achieve all the above features it was necessary to sacrifice some 
of the accuracy of the calculations. In the worst case, only three digits are reliable, 
although in most cases four and five are correct. Consequently, the results are as 
good or better than obtained with a slide rule and, therefore, in the majority of cases, 
quite satisfactory for engineering work. 
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Card Layout and Programming Sheet 

The programming sheet is a duplicate of the card layout except for the punch- 
ing field (fig. 1). When planning for the electronic calculator, it is possible to use 
the programming sheet directly. When completed, it will also serve as a chart for 
the key-punch operator. The various fields in the card are allocated as follows: 
Column 1-4 Card numbers indicating the sequence of a program deck. 
Column 5-7 Special codes or any other indicative information which may 

be included with the results. 
Column 8-11 Operating and transfer instructions which control the func- 
tions of the machine. 
Column 12-18 Reading instructions for the basic storage locations. 
Column 19 - 60 Fields of the seven five-digit storage locations, each of 
which includes a column to indicate a negative quantity. 
Column 61 - 80 Punching fields for four five-digit results. 

Storage Units 

Factor storage and general storage locations are indicated by the letters F 
and G, respectively. In order to maintain uniformity, all storage units can accept 
only five digits and a sign. The basic storage units which are involved in all mathe- 
matical operations are factor storage 1-3, the multiplier-quotient unit, and the coun- 
ter. For simplicity, the numbers contained in those locations will be denoted by 
F, M and C, respectively. 

In order to read a number into electronic storage, it is necessary to have the 
number punched in the appropriate field together with an X-punch in the correspond- 
ing read column (12-18). If the number is negative, an X is punched in the column 
denoted by a minus sign. The column is left blank for positive numbers. This ar- 
rangement is used to avoid double punching. 

The three general storage units are not involved in any calculations and can 
be used for storing factors and results. The contents of these locations are also 
available for punching. 

Storage units F2 and F4 are used for some calculations. When not in use, 
however, they are available as temporary storage locations but their contents can- 
not be punched. In particular, during addition, subtraction, multiplication, and div- 
ision both F2 and F4 can be used as intermediate storage locations. When square 
roots are evaluated, only F4 is available. Neither F2 nor F4 can be used for storage 
during the calculation of natural or inverse trigonometric functions, or during the 
table look-up operation. 
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Storage location Fl-3 is the basic unit in which the numbers to be operated 
on are contained. When the operation involves two numbers, as in addition, sub- 
traction, etc. , the other factor must be located in M. In all operations except div- 
ision, the number F, stored in Fl-3, is retained throughout the calculation. At the 
end of any operation the five -digit result is available in M and C. 

All storage positions retain their numbers until another number is read into 
them. The electronic counter, C, is the only exception; it is reset at the end of each 
card cycle. 

Operating Instructions 

Columns 8-11 control the operation of the machine. There are two main 
groups of instructions which can be given to the calculator. Columns 8 and 9 contain 
primarily the codes for arithmetic functions, whereas columns 10 and 11 indicate 
the transfer operations. In most cases the two functions are independent and can be 
performed simultaneously. 

The details of the various codes are given below: 

Arithmetic Functions 
Addition. Code 1 in Column 8. 

The operation performed is . 1(F+M). Both factors are assumed to be of 
order .xxxxx and the results of order x.xxxx, correct to four decimal 
places. 

Subtraction . Code 3 in column 8. 

The operation performed is . l(F-M). Both factors are assumed to be of 
order .xxxxx, and the result of order x.xxxx, correct to four decimal 
places. 

Multiplication. Code 2 in column 8. 

The operation performed is F* M. Both factors are assumed to be of 
order .xxxxx, and the result of order .xxxxx, correct to five decimal 
places. 

Division. Code 4 in column 8. 

The operation performed is . 1F/M. Both factors are assumed to be of 
order .xxxxx, and the result of order x.xxxx, correct to four decimal 
places. Should the result be greater than 9. 9999, a division overflow 
will occur. This will cause the machine to stop and the zero check light 
to go on. The card which caused the division overflow will be at the top 
of the stacker. 
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Square Root . Code 7 in column 8. 

The operation performed is JW. The number F and the result are both 
assumed to be of order x.xxxx. The maximum error can be ±. 0001. If 
it is known that g F S 1.0000, then 10 cards are required in order to 
evaluate the square root to within the accuracy stated above. The iter- 
ative process is based on the formula 



= J_( C + I, 
'n + 1 2 * n + c )» 



C n 

where F is the number whose root is required and C n is the nth approx- 
imation. In the above case the first approximation should be taken as 
0.1001, which has to be read into F2 from the first card of the square 
root sequence. This should be followed by nine cards, each punched with 
a 7 in column 8. If the order of magnitude of F cannot be specified as 
above, the first approximation in F2 should be 1.0000 and followed by 
13 iterations, making a total of 14 cards. At the end of the iteration the 
result is also available in F2, in addition to the usual locations M and C. 
Cosine . Code 6 in column 8. 

The operation performed is cos F. The angle F, in radians, is assumed 
to be positive, and of the order x.xxxx. The maximum value of F is ir . 
The result is of the order x.xxxx, and at least correct to three decimal 
places, except for angles larger than about 150°. The iterative process 
is based on the formula 

F 2 F 2 „ F 2 u 

cosF = l- I- {1-J_(1 - ±- •••)} 
1.2 3.4 5.6 

Eight cards are required to obtain the accuracy given above. The num- 
bers 01400 and 10000 should be read into M and F4, respectively, from 
the first card of the sequence. This should be followed by seven cards, 
each punched with a 6 in column 8. At the end of the iteration the result 
is also available in F4, in addition to the usual locations M and C. See 
Appendix 1 for calculation of other natural trigonometric functions and 
for calculations of cos F when F is negative. 
Inverse Sine . Code 5 in column 8. 

The operation performed is 1/F sin F. The factor is assumed posi- 
tive, of the order .xxxxx, and is limited by the conditions that SFSl/y/27 
The results are of the order x.xxxx and the maximum error is .0001. 
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The iterative process is based on the formula 

-Isin-^l-H 1 2 F 2 {1+JLf 2 (1+ 5 2 F 2 .. )} 
F 2.3 4 ' 5 6.7 

Eight cards are required to obtain the accuracy given above. The num- 
bers 01500, 00015, and 10000 should be read into M, F2, and F4, res- 
pectively, from the first card of the sequence. This should be followed 
by seven cards, each punched with a 5 in column 8. At the end of the 
iteration the result is also available in F4, in addition to the usual loca- 
tions M and C. See Appendix 2 for calculation of other trigonometric 
functions and for calculation of sin ~ F/F when l/y/Ts F § 1 or -1 ^ F = 0. 



Transfers 



It is possible to transfer a factor from any storage unit into Fl-3 before the 
calculations take place. Consequently, the transfer and the operation on that number 
can, in general, be performed during the same card cycle. The two exceptions are 
M and G4, which should not be used in calculations on the same card cycle as the 
transfer takes place. After each calculation, i.e. , during the punching cycle, the 
results can either be punched or transferred from the counter into any one of the 
storage locations. Since at that time M and C contain the same number, the transfer 
between these two units is not required. It is also not possible to transfer from C 
to Fl-3, but the transfer from M could be used instead. The numbers transferred 
from C are available on the following card cycle. 

The transfer codes are given in detail below: 

Transfer from F2 to Fl-3. Code 4 in column 10. 

Transfer from F4 to Fl-3. Code 5 in column 10. 

Transfer from Gl-3 to Fl-3 . Code 1 in column 10. 

Transfer from G2 to Fl-3. Code 2 in column 10. 

Transfer from G4 to Fl-3. Code 3 in column 10. 

The number transferred can be used for calculations on the same card cycle 

only if the shift code (code 9 in column 9) is not punched in the same card. 

If the shift code is present, the number in Fl-3 will be shifted one position 

to the left before the arithmetic operation takes place. This feature may be 

of advantage in some cases. 
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Transfer from M to Fl-3. Code 3 in column 10 together with code in col- 
umn 8. The remarks of the previous paragraph apply here as well. In ad- 
dition, since this transfer requires punching in the operation column 8, 
double punching will be required in order to perform an arithmetic opera- 
tion during the same card cycle. 
Transfer from C to F2. Code 9 in column 11. 
Transfer from C to F4. Code Y (12th position) in column 11. 
Transfer from C to Gl-3 . Code 6 in column 11. 
Transfer from C to G2 . Code 7 in column 11. 
Transfer from C to G4. Code 8 in column 11. 

Transfer to Fl-3 takes place at the beginning of the calculating time, while 
transfer from C occurs during punching time. Consequently, both codes can be used 
independently on the same card. If a code in column 11 is used, the transfer will 
take place when the card containing the code reaches the punching station. 

Shift. Code 9 in column 9. 

During addition, subtraction, and division the decimal point is shifted one 
place to the right. This is necessary, since the order of magnitude increases dur- 
ing those operations. When the numbers which are involved in the calculation are 
small, it may be desirable to maintain the decimal point at the same place. This 
can be done by shifting the result one position to the left. If a shift is used during 
multiplication, the order of magnitude is reduced by one. 

The above can be summarized symbolically as follows: 

Add and shift . xxxxx + . xxxxx = . xxxxx 

Subtract and shift . xxxxx - . xxxxx = . xxxxx 

Multiply and shift . xxxxx x . xxxxx = . Oxxxxx 

Divide and shift . xxxxx_i_ . xxxxx = . xxxxx 

In the case of addition and subtraction the results are correct to five decimal 
places. Only five digits are available for the result of multiplication, and the first 
decimal place shown as zero does not appear in the result, which is correct to six 
decimal places. In the division, the last digit is not half adjusted, i.e., the maxi- 
mum error is .00001. The result cannot be shifted during the evaluation of square 
roots, cosines, or inverse sines. The use of the shift code in that case will be in- 
effective. In some cases shift can be applied together with a transfer operation. 
This feature is normally used for transfer of numbers from M to Fl-3, although the 
shift can also be made during transfer from G4 to Fl-3. The above is equivalent to 
direct multiplication by ten. 
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Test for Sign of Arguments . Code Y (12th position) in column 8. Code 8 in column 9. 

The calculation of cosine and inverse sine can only be made for positive argu- 
ments. It is therefore essential to be able to determine the sign of the argument be- 
fore the evaluation of the function is started, and to change the number to its absolute 
value. 

If code Y in column 8 is then used on the following card cycle, the factor con- 
tained in M will be changed to its absolute value. Moreover, if the original value of 
M was negative, one of the pilot selectors will be picked up and remain transferred 
until a special code is used to bring it back to normal. The latter is code 8 in column 
9, which can be used on any subsequent card. This code also enables the subtraction 
F-M to be performed on the following card cycle if the original argument in M was 
negative; otherwise no operation takes place. 

To illustrate the above, consider the calculation of cos "A. As a preliminary 
step the argument A is transferred to M and tested for its sign by means of code Y in 
column 8. On the next card cycle the value |A| will be given in M. A series of cards 
will now follow in order to perform the calculation which is only possible for positive 
arguments (see Appendix 2). At the end of this process the result, say X, is again 
inserted in M, and the value ir is read into Fl-3. The code 8 in column 9 is then 
given. If originally A was negative, subtraction will be performed on the next card 
cycle, giving {it - X) as a result; if A was positive, no operation will take place and . 
X will remain as the result. When code 8 in column 9 is used, it is assumed that both 
numbers M and F are of order .xxxxx and the result is also of the same order, cor- 
rect to five decimal places. 

The test for negative arguments can also be made before the extraction of a 
square root. In this case a test division could be programmed in such a way as to 
cause overflow if the factor was negative. As a result the machine will stop, indica- 
ting that the square root operation was invalid. 

The only limitation for the above test is that no result may be punched on any 
intermediate card until the particular calculation is completed. In other words, there 
should be no card containing code X in column 11 (punching code) between the card 
having code Y in column 8 and the card with code 8 in column 9. It is unlikely that 
such a case would be required in practice. 
Punching. Code X in column 11. 

All punching is done on blank cards. Consequently, the code which actuates the 
operation should be given on the preceding card. When code X in column 11 is given, 
the contents of Gl-3, G2, G4, and C are punched into their respective fields on the 
following card. 

Sometimes it may be necessary to include identifying codes together with the 
punched results. This may be done by gang punching into the special codes from the 
preceding or succeeding card, depending whether the gang punching is wired from 
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second or first reading. If several results are interspersed in the same deck, it may 
be convenient to gang punch the card number of one of the adjacent cards. During the 
punching operation factors cannot be transferred from C into any of the storage loca- 
tions. 

Table Look-Up. Code X and 3 in column 8 together with codes 8 and 9 in column 9. 
Code X in column 8 together with code 8 in column 9. 

The table look-up has two distinct types of cards. The first card of each table 
prepares the machine for the table look-up operation. This card has codes X and 3 
in column 8 together with codes 8 and 9 in column 9. The remaining cards contain the 
actual table and have code X in column 8 together with code 8 in column 9. Before the 
start of the operation the searching argument is placed in Fl-3. The layout of the 
table cards is as follows: the arguments are in field M, the function in field F2, and 
the first order differences in field F4. During the operation, the argument in M is sub- 
tracted from the searching argument in Fl-3 and the fields F2 and F4 are read in. If 
the difference between the arguments is negative, the operation is repeated for the fol- 
lowing card and the process goes on until a zero or positive difference is obtained. 
Subsequently, the remaining cards of the table pass through without being read and sub- 
traction is suppressed. Consequently, at the end of the table look-up the searching ar- 
gument is in Fl-3; storage unit M contains the difference between the searching argu- 
ment and the first argument of the table which is smaller or equal to it; F2 and F4 con- 
tain, respectively, the function and first order difference for the latter argument. 

The table must always be arranged in a descending order of arguments. No 
reading instructions (columns 12-18) should be included in the table cards. In order to 
prepare the machine for the table look-up, it is necessary to perform a subtraction 
which gives a negative result. This is obtained by reading the number 99999 into M. 
Consequently, the first table card should consist of the appropriate codes (X-3 and 
8-9)99999 in field M, reading instructions for M (column 12), and possibly the search- 
ing arguments and reading instructions in field Fl-3. The subtraction during the table 
look-up is of the type F-M or .xxxxx - . xxxxx = .xxxxx, correct to five decimal places. 
The table look-up is possible for positive arguments only. 
Wiring Diagrams 

The wiring diagrams are shown in various sections for ease of reading. Fig- 
ures 2, 3, and 4 show the wiring of different parts of the punch panel. In a few cases 
the same split-wire is given on more than one diagram. The connections for the cal- 
culator selectors on the calculate control panel are given in symbolic form in figure 
5. The meaning of the symbols is as follows: 

(23J A number enclosed by a circle - suppression hub of a particular program 
step. 
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(g) A letter enclosed by a circle - suppression hubs of a suppression group. 
\3l) A number enclosed by a rectangle - one program outlet of a particular pro- 
gram step. 
S Suppress without balance test. 

SUP + Suppress on plus balance. 
SUP - Suppress on minus balance. 
BTSS Balance test for step suppression. 
BTSPU Balance test for selector pickup. 
RI, RO followed by the symbol of a storage unit denotes read into, and read out of 

the indicated storage unit. 
RI + Read into counter + 
RI - Read into counter - 
M + Multiply + 
M - Multiply - 
DIV Divide 
Z.CH. Zero check 
1/2 Half adjust 

RUI, RUO followed by a number denotes read unit into, and read unit out of the posi- 
tion indicated. 
ET followed by a number denotes number obtained from emitter control. 
RO Read out of counter 
R&R Read out and reset counter 

The individual program steps which are a guide to the wiring of the calculate 
panel are given in tabular form in figures 6(a) and 6(b). Since depending on the oper- 
ation codes, etc. , different operations may be performed on the same program step, 
all the possible combinations are included in the table. Each line indicates the codes 
corresponding to the particular operation steps and the calculator selectors which are 
transferred at that time. The program steps are combined into a number of suppres- 
sion groups, which are referred to in figure 4. Suppression group B is different from 
all others and is indicated by *. The three program steps included in that group are 
always operative, except when calculator selector 3 is transferred, i.e. , either when 
code 3 or code 5 are used in column 8. 

The wiring of the punch panel presented considerable difficulties, and it was 
not possible to avoid all the back circuits owing to a shortage of selectors. Consequent- 
ly, three small rectifiers had to be added and these are indicated by the standard sym- 
bol. Two types of rectifiers were tried and both gave satisfactory results: the selenium 
type, consisting of three plates, each one inch square; and germanium diodes, only 
about 3/8 ''long and 3/16" in diameter. 

Another effect of the shortage of pilot selectors is some unusual features used 
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in the wiring. In particular, both double punch detection and product overflow are 
used as delayed pickups. 

The programming is such that there is only just enough time to perform all the 
calculations during one card cycle and the unfinished program is wired to stop the 
machine as a safety precaution. 

Appendix 1 

The calculator is capable only of evaluating cosines of positive angles. If 
other trigonometric functions are required, they can be obtained from the following 

formulae: . — 

sin X = /1-cos X 
sin X 



tan X = 



cosX 



In order to ensure that the iteration is performed for positive arguments only, 
the absolute value of the argument should be obtained as a first step. The required 
function is then calculated using the above formulae. At the end of the calculation, one 
of the following subtractions is performed if the factors were initially negative: 

- sin X 

- tan X 

2 cos X- cos X 

Appendix 2 

A large number of iterations is required to obtain 1/X sin" X accurately for X 
close to unity. In order to reduce the calculating time, the value of X has been limi- 
ted to 1/ v2, corresponding to an angle of 7r/4 or 45°. 

Let Y be any positive argument such that 

0IYS1 

and let sin " 1 Y = A=-|— +B 

If B is now considered as the unknown quantity it follows that 





4 " B --4- 


Let 


B = sin - 1 X 


hen 


- 1 *X i 1 



V2 /2 
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A = sin Y 


for ^ Y s 1 


X = * (Y - VT^Y 2 ) 
/2 


B = sin _1 X 




A=-J-+B 




A = cos " Y 


for s Y i 1 


x = JL (Y-VUY 2 
/2 


~) 


B = sin _1 X 




A=_f-- B 




A = tan _1 Y 


for s Y i 9. 8 


x _ 1 X-l 




/2 ^+ 1 




B = sin _1 X 





The relationship between X and Y is as follows for the different inverse trigono- 
metric functions: 
To obtain 

find 

giving 
and 
To obtain 
find 

giving 
and 

To obtain 
find 

giving 

and A = -J— + B 

Due to the limited calculator capacity, the maximum value of tan " Y has to be 
limited to just over 84°. 

To ensure that positive arguments only are used, the procedure is the same as 

that given in Appendix 1 , except that one of the following subtractions is performed at 

the end: 1 

- sin Y 

IT - cos Y 
ir ■- tan _1 Y 
The above subtractions can be modified if different quadrants are required. 

Appendix 3 

For most calculations it is arbitrarily assumed that the numbers are less than 
unity and the various powers of 10 are obtained by shifts to the left or right. In the 
case of some operations, however, it is important to locate the decimal point at the 
place determined by the programming. It is therefore suggested that the following 
scheme be used throughout. 

. xxxxx + . xxxxx = x. xxxx without shift 

= . xxxxx with shift 

. xxxxx - . xxxxx = x. xxxx without shift 
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= . xxxxx 

. xxxxx X . xxxxx = . xxxxx 

= .Oxxxxx 
. xxxxx -4- . xxxxx = X. xxxx 

= .xxxxx 

Vx.xxxx = X. xxxx 

cos (x.xxxx) = X. xxxx 

1 . -1, 

sin 



(.xxxxx) 
References 



' (. xxxxx) = x. xxxx 



with shift 
without shift 
with shift 
without shift 
with shift 
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NUMERICAL INTEGRATION OF HYPERBOLIC PARTIAL DIFFERENTIAL 
EQUATIONS WITH CONSTANT COEFFICIENTS AND WITH CHARACTERISTIC 
INITIAL CONDITIONS ON THE 604 ELECTRONIC CALCULATING PUNCH 

N. Arne Lindberger 

International Business Machines Svenska AB 
Stockholm 

This paper is concerned with the numerical solution by means of punch cards 
of the equation 

9 ^ +a -i£- + b-!4- + c0 = O, 
3£ 817 3 £ 3 77 

where a, b, and c are constants. The initial values of are given along two orthog- 
onal characteristics. 

In section 1.1, by transformations in the dependent and the independent vari- 
ables, the equation is reduced to the form 

3x 3y 

In section 2.1 by approximation to a difference equation, a step-by-step pro- 
cedure is derived, whose error and stability properties are analyzed in sections 2.2 
and 3.1. On the basis of these results, section 4.1 describes a modification of the 
step-ahead process which has been set up for automatic integration on the 604 for the 
particular starting conditions z(x,0) = 0; z(0,y) = e y -l. The corresponding values of 
have also been calculated for the case of a=b=l, c=0. 

The machine set-up is described in sections 6 and 7. 

In sections 5.1 and 5.2 the general solution of the problem is given as two 
different series expansions which have been used for an accuracy check on the results 
of the 604 runs. 

1.1 By means of a linear transformation in the independent variables, any linear 
homogenous hyperbolic differential equation with constant coefficients can be reduced 
to the canonical form 

3 2 +a ±0 + h 9* 



+ a-ii_+b U 4- C = 0, (1.11) 



3|3n 3| 877 

where a, b, and c are constants. 

The boundary conditions are assumed to be given along two characteristic 
curves, belonging to different sets; that is, £ and 77 equal constants. 

The partial derivatives of the first order can be removed from equation (1.11) 
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by the transformation 



0(«,ij) = z(S,i? )'e 



• ~ - *>£ - ar? 



(1.12) 



Equation (1.11) then reduces to 



a 2 z 



= (ab-c)z 



3£ 877 

The trivial case of c=ab will be disregarded. 
At last 



(1.13) 



x = a g c • £ ; y= K77 , 



(1.14) 



where Kis an arbitrary constant, will transform (1.13) into 



8 2 z 

3x 3y 



= z 



(1.15) 



which is the form that is going to be treated numerically. 



1.2 Formal integration of (1.15) over the domain indicated in Figure 1 yields 3 



xy 



z M = z p + z Q - z Q + / / z(a, /3) da d£ , 



00 



(1.21) 



where a, /3 are integration variables. 

y. 



M(x,y) 




Figure l 

2.1 For an approximation to finite variables the first quadrant of the xy-plane is 

covered with a net of equal rectangles by the straight lines 

x= n ax, n = 0, 1, 2, 3, ... 
y = mAy, m= 0, 1, 2, 3, ... 
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m + 1 



m 



y 

n, m+1 




n+1 , m+1 


k=Ay 




P M 
Q 


U 1 

n+l,m 




u 
n, m 










h=Ax 





-> x 



n 



n+1 



Figure 2 



The equation (1.21) is now applied on a function u n defined at the nodes of the mesh 
and intended to be an approximation to z(x,y). The integral in (1.21) is approximated 
by 

(2.11) 



hk 



UK/U, .. + U ., +U ,+U v 

-£— ( n+1, m+1 n+l,m n, m+1 n,m) , 



where 



h = ax, k=A y, u n m ~ z(x,y) , 



i.e. , the integrand is treated as a constant which is equal to the mean of the nodal 
values. 

Insertion into equation (1.21) now yields 



u n+l,m+l = < u n+l,m + u n,m+l> < 1+ 6 > " u n, 



m 



6- hk/2 
l-hk/4 



(2.12) 
(2.13) 



Incidentally, this formula is not only a first approximation: it also represents the 

4 
limes of a Picard's iteration process (see section 5.1) on equation (1.21) with the 

integral approximation (2. 11). 

Equation (2.12) is the equivalent of (1.15) in the finite difference domain. Of 
course, other approximations could have been obtained by other means of approxima- 
ting the integral. Before applying equation (2. 12) to numerical integration, its error 
properties will be investigated from various aspects. 



2.2 The function u n m defined by (2.12) and (2.13) is supposed to be an approxi- 

mation of the function z(x,y) defined by (1.15). When inserting z into (2.12), it will 
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be necessary to add an error quantity t(x+h, y+k) to the left-hand side 

z n + l, m + l + Wl,m + 1 = ( Vl, m + \ m + l> (1+6) " \ m (2 ' 21) 

By using the Taylor's series for 

z n + l , m H z(x+h ' y) and Z n, m+l H z(x ' y+k) 
and expanding 6 of (2. 13) in a geometrical series, the right-hand side of (2. 21) will be 

,, hk h 2 k 2 . , . 3z A _ h 2 3 2 z ^ '[ 

( i+ t- + -5— + --- )(z + h 9x" + nr Tl? + i 

h 3 3 3 z h 4 3 4 z ^ 17 . 3z k 2 3 2 z 
+ — — — + — — + . . . . + z + k + — + 

13 3x 3 (4. 3x 4 3y \2, dy 

+JL--g!g. + JL. 9 z +...) -z (2.22) 

io ^a \a - 4 ; n,m n,m 

H 3y 11 3y ' 

When multiplying with hk/2, according to its definition in equation (1. 15), z shall be 

o 
regarded as 3 z/3x3y. The same is valid for the partial derivatives of z. When 

2 2 4/22 

multiplying with h k /8, z is set equal to 3 z/3 x 3 y. 

With these rules, the expansion of (2.22) gives 

z n,m + h N m * k \m + * { h 2 9 \m + 2 tt 9 2 n>m + 
2 3 ? 

«3Z 10 3Z O 3 t; 

+ k 2 "> m } +4- { h 3 ,. n r m + 3 h 2 k Z n,m 



3 o 4 
9 3 z <i g° z i 4 3 z 

+ 3 hk 2 n i m + k d d Z n,m + -|_ h 4 n,m + 

3x3y 2 3y 3 U X 3x 4 

4 4 4 4 

+ 4 h 3 k a \m + 6h 2 k 2 9 V m + 4 hk 3 9 Vj + k 4 % 

3x 3 3y " 3x Z 3y 2 3x 3y 3 Ty 3 

+ hk (h 2 9Vm_ 2 jVm ) , x 

3x 3y 3x 3y° 

This is nothing but the Taylor series expansion of the left-hand side z^ j m+1 = 
z(x+ax, y+Ay) plus an error term of the fourth order. Thus it has been shown that 
when using the recurrence formula (2.12), (2.13) for computing z, defined by (1.15), 
in each step of the process a truncation error 

t ( X ,y) = ^dLAl { ( Ax) 2 _J^_ + ( Ay) 2 _A } (224) 

12 3x 2 3y^ 
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is committed. 



3. 1 After establishing the degree of approximation to the partial differential equa- 

tion when using the difference equation (2.12), its stability properties will be inves- 

5 
tigated. Using a method attributed to von Neumann , assume a solution in separated 

variables 



u 



= ^. 



im0 



(3.11) 



Insertion into (2.12) yields 



= 



e 1 * (1+5 )-l 



(3. 12) 



I )3| always equals 1 independently of and 6. This means that the choice of step 
lengths A x and A y will not influence the stability of (2. 12). There will be no serious 
growth of errors, nor will they be reduced. 

This simple criterion is confirmed by an €- diagram . Assuming an isolated 
error e to have occurred in the computation of u the growth of this elementary 

quantity is illustrated by Figure 3. 



m+M 





e(l+M6) 


e[l+(3M+l) 6 ] 


e[l+(5M+2) 6 ] 




e[l+(2MN+M+N) 6 ] 
















m+2 





e(l+2 6 ) 


e(l+7 6 ) 


e(l+12 6 ) 




e[l+(5N+2) 6 ] 


m+1 





e(l+6) 


e(l+4 6 ) 


e(l+7 6 ) 




e[l+(3N+l) 6 ] 


m 





6 


e(l+6 ) 


e(l+2 6) 




e(l+N 6) 

























n 


n+1 


n+2 




n+N 



Figure 3. e-diagram for the step-ahead process (2.12). 

Disregarding all powers of 6, after N steps in x and M steps in y, a single 
error will have grown to e[l+(2MN+M+N) 6 ]. As M, N «CMN and 6 ^ AX 2 Ay , 
the expression for the accumulated error reduces to e(l+X* Y) where X = Na x, 
Y= MAy. 

Suppose that an error e has occurred in the first step of the first column. In- 
tegrating over, say, a square with side length 3 units in the first quadrant, when the 
upper right corner (3, 3) of the square is reached, the accumulated error due to the 
single disturbance will be lOe. Although not prohibitive, this rate of error accumula- 
tion is undesirable. Without increasing the truncation error (2.24), the situation can 
be remedied by using another method of step-ahead integration. 
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4„1 An analysis of the step-ahead process along the straight line x = (n+1) &x 

will show that it is equivalent to the following summation formula: 

r ^ / u • + u 1 • 1 \ U n+1 ,0" u n,Oi s (a ti\ 

Vl,m + l = u n,m + l + [ .^ < n,x n+l f i-l) + ^ *- ] 6 (4.11) 

Here for the first time, the initial values along the characteristic y=0 are entered in- 
to the equations by the last term inside the brackets in (4.11). 
With the definitions 



m 



u_ . + u_ 



*n+l,m = £< n > 1 n+1 ' i ' 1 ^ + a n + l,0 



(4.12) 



a. 



U n+1,0 - u n,0 



'n+l,0~ g 

a summation step from m to m+1 can be broken up into the following operations: 



v n+l,m+l =u n,m+l + U n+l,m 



a n+l , m+1 a n+l , m + v n+l , m+1 



u n+l , m+1 u n, m+1 + * ^n+l , m+1 



(4.14) 
(4.15) 
(4.16) 



The growth of an isolated error in u during this process is shown in Fig- 

n n, m 

ure 4, where powers higher than 6 have been disregarded. 



m+M 





e6[l+(M-l)6] 


e 6[3+(5M-3) 6 ] 


efi [5+(13M-5)6] 




€ 6{2N+1+[M(2N 2 +2N+1)-(2N+1)]6 } 
















m+2 





e6(l+6) 


e 6(3+7 6) 


e (5+21 6 ) 




c 6 [2N+1+(4N Z + 2N+1) 6 ] 


m+1 





e 6 


e 6 (3+2 6 ) 


e (5+8 6 ) 




e 6 (2N+1+2N 2 6 ) 


m 





e 


e(l+6) 


e (1+6 f 




€ (1+6) N 

























n 


n+1 


n+2 




n+N 



Figure 4. e -diagram for u derived from (4.14)-(4.16). 

Here an error propagates mainly along a horizontal line in the diagram. As N5^C 1, 
it will practically not grow. Above the main line of disturbance, the error will pro- 
pagate according to the expression 

e 6 {2N+1+[M (2N 2 +2N+1) - (2N+1)] 6 } 
With the same approximations as in section 3.1, this expression reduces to 
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e 6 • N • (XY + 2) 

For X=Y=3, N=300 this gives 0. 165e, which is a considerable improvement compared 

to the previous process in section 3. 1. 

The corresponding e -diagram for a _ shows the same trends, but in limes 

n,m ' 

values are doubled as per the diagram in Figure 3. This is of no consequence, how- 
ever, as the function wanted is u , not a 

' n,m' n,m 

It should be stressed that the developments of section 4. 1 are equally valid 
if n and m are permuted, i.e. , the summations are performed in parallel to the x- 
axis. The choice is a matter of convenience with respect to the initial functions. 

4.2 Integration of equation (1.15) with respect to y for a constant x yields 

9Z a (X ' y) = 7 z(x, 0) d/3 + X '(x) (4. 21) 

3x o 

X (x) is the initial value of z along the x-axis. 

If in (4. 11) 6 is approximated by Ax A y/2 and the equation is slightly re- 
written, an equivalent of (4. 21) in the finite domain will appear 

u n + l,m+l- u n,m + l = Ay [ m ^ u + %+1 } + Vl,0 " \o ( 4 .22) 
ax i=l Ax 

The three main terms of (4.22) correspond to the three terms of (4.21). 

5. 1 Under the usual assumptions of continuity and differentiability the higher de- 

rivatives of z with respect to x can be written as 

iZ^ = j a^ffl) dg + x (p) (x) (B.U) 

3x p 3x p_1 

Thus every new derivative can be calculated from the preceding one. After introdu- 
cing the initial values along the y-axis, z(0,y) = i//(y), where xf/(0) = x (0), the solu- 
tion of (1.15) with the specified initial conditions can be expanded in a Taylor series. 

Z (x,y) = z(0,y) + x **«>>* + A 1*&& + A £*&& (5 .12) 

3x £ dx' y. dx 6 



After insertion of (5.11) and observing that 

• V P 

P=l |£ 



s JL X (P) (0) = x (x) -x(0), (5.13) 
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(5. 12) yields B Q 

z(x,y)=tf/(y) + X(x)-X(0)+ S -£-{ / /•••/ ^0 X ) d^ dj3 2 . . . dfi + 

p=l li- • . 

q=l Li 

y 

For p=l the multiple integral is to be interpreted as the single integral / VOi) d/3 1 
p-1 

and the sum 2 a s zero. 
q=l 

Within the radius of convergence of the Taylor series, equation (5.14) repre- 
sents another means of numerical solution of equation (1.15) with the characteristic 
initial values x(x) and ty (y) along the axes. As x and y are permutable in (1.15), per- 
mutation of x and y, x and \j/ in (5. 14) will give another resolution formula for z(x, y). 

5.2 A series expansion without the restrictions of differentiability on the deriva- 

tives of one of the initial functions can be obtained from (1.21), using Picard's sub- 
stitution method . Assume a first approximation: 

z (x,y)=»My) + x(x)-x(0) (5.21) 

Insertion into (1. 21) gives a new approximation 

y x 

z.(x,y) = z n + x JiK/3,) ft + y / x(«i) da - xy • X (0) (5.22) 

i u i i 1 1 

Repetition of this process will give a converging series 

z(x,y)=^(y) + x(x)-x(0)+ 2 -^- • / / .. . / #(%) d0, d/3 ... da + 

P =iL£ooo li p 

°° p x a p a 2 °° p p 

p=iL£oo o^ 112 p p=i ( \pr 

This expansion represents the unique solution of (1.15) with the specified initial con- 
ditions . The only restriction on x(x) andj//(y) is that they shall have a first deriva- 
tive throughout the region of integration. 

If x(x) is supposed to be expansionable in a Taylor series, the identity be- 
tween (5.23) and (5.14) is easily established. Both these formulas have limited numer- 
ical importance because in practice the multiple integrals can seldom be explicitly per- 
formed. 

6. 1 For a test setup on the 604 a specific set of constants and initial values have 

been chosen. In equations (1.11), (1.12), and (1.13) a and b have been set equal to 1 
and c equal to zero. 
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Initial values: For x = 0, 


0(O,y) = l-e" y 


For y = 0, 


0(x,O) = O 


By the transformation 


= z e" (x+y) 


(1.15) will reappear 


/*~ 



Therefore, the setup problem is the following*: 

_i!l_ + ii_ + M_ = o (6.ii) 

9x Sy 3x 3y 

(6.12) 

(6.13) 

3x 3y 
with the initial values 

Forx=0, z(0,y) = ^/(y)= e y -l (6.15) 

Fory=0, z(x,0) ■ x (x) = 

Insertion of ^and x into either of equations (5.14) or (5.23) gives the solution of the 

problem in the form of the series expansion 

2 2 
z(x, y) = e y -l + (e y -l-y) x + (e y -l-y-j|-) J- + . . . (6. 16) 

After rearranging the order of the terms, this expression can be written 

z(z,y)=e x+y - S j£ s ■£- (6.17) 

p=0 LE q=0 \1 

and in this form coincides with one of the mentioned Bessel expansions. From 
(6. 13) and (6. 16) or (6. 17) it is seen that 

OS0<1 (6.18) 

This limitation on is going to be used as a check in each computation step. 
Equation (6. 16) is used as a manual check when the machine integration is completed. 

6.2 The main computation formula is (4.11), where the term u - n -u n /» 

disappears owing to (6.15). Every step in the calculation is checked against equation 
(2.12). 

One run on the 604 represents the integration along a straight line, say, 
x= (n+1) ax in the xy- net, Figure 2. This line will be called column n+1. According 

* This case has been treated by C.G. Allander in his doctor's thesis "Untersuchung 

des Adsorptionsvorganges in Adsorbentenschichten mit Linearer Adsorptionsiso- 

therme," Royal Institute of Technology, Stockholm, 1953. Using the classical Rie- 

mann integration method for hyperbolic equations, Allander expresses the solution 

of equations (6.11), (6.12) as y where I n is the Bessel 

= e~ x / e _t I ft (2 Vxt) dt, u 


function of order zero with pure imaginary argument. After expanding the integrand 

in a series and integrating, he arrives at the expression 

0=1 _ e -(x + y) s(j)ll ffi <2V*y) 
which is identical with equation (6.17) after multiplication with e x+y . 
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to (2.12) and (4.11) the values u - - are entirely determined by the corresponding 
values u - one column back, together with the first value u - Q = in the new 

column, ax and Ay have both been chosen equal to 0.01. This makes 5 in (2. 13) equal 

-4 
to 0. 50001 • 10 . This is fed from the digit emitter. 

The exponential function is calculated for each integration step with the equa- 
tion 

e -< x+ y>n + l,m + l = e -Ay_ e -C**y) n+1>m , (6.21) 

where e" Ay = 0. 9900498. It is fed from the digit emitter as 0. 99005. 

e~^ x+y Wl,m is available from the previous step. Only the first value of 

a new column e~^ x+y ^n+l, has to be read from the card which is used for starting the 
integration process. 

7.1 In the test runs the integration has been carried out to 3 units in y. This 

means that there are 300 cards for the values u n . Before each run these master 

cards are interspersed with blank cards, which are going to be punched in columns 

18-25 with the results u - of the actual run and then become the master cards of 

n+l,m 

the next integration run n+2. 

The blanks are also punched in columns 27 - 33 with ^ ^which is the 
ultimate information wanted, and in columns 34 - 41 with e _ ( x+ y) n+ i m > which is pun- 
ched out only to be read into the machine again from the second reading station and 
used according to (6. 21). In this respect the card is used as a moving auxiliary stor- 
age, to remedy the limited storage capacity of the 604. 

The mentioned quantities are punched from the counter and general storage 
of the 604 as specified in the flow diagram, Figure 5, below the heading "1st Card 
Cycle." The blank card is also punched with the number m for identification in col- 
umns 5-7 and with u in columns 9-16. Both these quantities are read from the 

n, m 

preceding master card m in the second reading station, u is also to be read into 
the machine again during the second card cycle when the former blank card reaches 
the second reading station for use in the check formula (2.12). The electronic stor- 
age does not suffice for keeping this quantity over the two machine cycles which to- 
gether constitute one computing step. 

The only original information in the blank cards is the x-coordinate number 
n+1 which is gang-punched into columns 1-3 before the merging with master cards. 
The master cards are not punched during the actual run, only during the preceding 
run when they were blanks. They control various machine functions with an X-punch 
in column 80. The blank cards have to be punched with this X in order to be distin- 
guished as masters in the next integration run. 
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7. 2 The processing of information in position and time is illustrated in Figure 

5, Flow Diagram of Type 604 Operations, and in Figure 6, Card Diagram of Type 604 

Operations. 

The punching of a blank card m , described abov°, occurs in the beginning 
of a new computing step. The results of the previous computing step are punched out 
simultaneously with the reading of new data from master card m+1 for use in the first 
of the two calculating cycles in the new computing step. u n ^ is read through punch 
selectors 1,2, transferred by the X-punch in column 80 of the master card, into fac- 
tor storage 3,4, assigned together. The decimal point convention is such that u n m+1 
will have 5 decimals and 3 integer digits. 

The constant 6 = 50001. 10 is sent from the digit emitter to the multiplier 
quotient unit through punch selector 3, which is also controlled by the X-punch in col- 
umn 80. 

Selectors must be used for reading into various storages because in the sec- 
ond card cycle other data have to be read from other sources into the same units. 

The symbolic Memory Contents 1 in the flow diagram, Figure 5, indicates 

that in addition to u - other information is stored since the previous calculating 

n, m+1 

step, u 1 is stored in general storage 3,4 with the same decimal convention as 

for u H ! The quantity a . , defined in (4. 12), is stored with 3 decimals and 5 

n,m+l ^ J n+l,m' 

integers in factor storage 1,2. 

7.3 The operations during the first calculating cycle are shown in the flow dia- 

gram. The calculation is performed according to formulas (4.14) - (4.16), leaving the 
new values of a 1 * and u 1 - in the same storages as the old values, shown 
by Memory Contents 2 in the flow diagram. In order to economize with storage, the 
check formula (2. 12) is rearranged to 

u =v < , (1+6) - u . , , (7.31) 

n,m n+l,m+l v ' n+l,m+l ' ' 

where the bar over u indicates that it represents a back solution from the newly 

formed values of u n+1 | m+1 and v n+1>m+r 

u is left in the counter after the first calculating cycle. The operations 

n, m 

described are performed on program steps 2-22. Unlike the remaining program steps 
they are not suppressed, because calculator selector 2 is picked up by the X-punch in 
column 80 „ 

During the second card cycle the master card m+1 is moved to the punching 
station. No punching occurs, as the punch suppression was activated during the pre- 
vious card cycle. A new blank card m+1 has entered the first reading station without 
being read. The previous blank card m , which is now carrying information, is read 
from the second reading station, as shown in Figure 5. This results in the transfer 



62 



n+1 m 



Starting value for new column 

-(x+y) r 
e 



card 



; n+l,0 



0000 

I 2 3 4 

1111 

2222 

1 2 2 4 

3333 



0000 

5 6 7 3 

1111 

2222 

5 G 7 5 

3333 



0000 

9 10 II 12 
1111 

2222 

9 10 II 12 

33 33 



0000 

13 14 IS 16 
1111 

222 2 

) 1) 15 11 

3333 



0000 

17 13 19 23 
1111 

2222 

17 IS 19 2D 

3333 



0000 

21 22 23 24 
1111 

22 22 

71 22 23 24 

3333 



000(1 

23 28 27 26 
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2222 

25 26 27 26 
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0000 
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2222 
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3 333 



0000 

45 45 47 48 
1111 

2222 

45 46 47 48 
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00 

65 66 67 68 
1111 



222 2 

65 66 67 -i 

33 33 



0000 

59 70 11 72 
1111 



2222 

69 '0 71 72 

3333 



run 

00 

73 74 75 75 
1111 

2222 

73 74 75 76 

3 33 3 



0000 

77 79 79 60 
1111 



2222 

77 78 79 60 

33 33 



n+l^m 



000 

12 3 4 

1111 

2222 

12 3 4 

333 3 



0000 

5 6 7 8 
1111 

222 2 

6 7 8 

3333 



1111 




and CR card m 
punched from D. E. « 



2222 



2222 



0000 

53 54 55 56 
1111 

2222 



0000 

57 5» 59 66 
1111 

2222 



0000 

81 82 63 64 
1111 

2222 



Numbers in brackets not used when calculating column n+1 

-(x+y). 



0000 

55 66 67 68 
1111 

2222 



,0 000 

69 70 71 72 
1111 

2 222 

69 70 71 72 

3 333 



0000 

73 74 15 75 
1111 

2222 

U /4 75 .5 

3333 



0000 

77 76 79 80 
1111 

2222 

71 78 79 63j 

3 3 3 3 



i 



card m+1 



1 /U^^u 1/ ii . "^ x+y 'n.nv-L All information except n 
n_ m+1 ( n-l,m+l) n, m+1 ( y n,m+l e ^ **) printed on previous 804P 




'run 



0000 

33 34 35 36 
1111 



'Read in:o 



2222 



0000 

37 38 39 40 
1111 



2 2 2 2 



mm^m 



0000 

41 42 43 44 

1111 



■*2 rt-< 



00000000000000000000 

45 45 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 



1111 






tin 



mi 

53 54 55 56 

3333 



1111 



2 

57 53 59 60 

333 3 



oooooo o old ooo a o u a 

65 66 67 69 63 70 71 72 73 74 75 76 77 79 79 80 



1111 



&m$nm 



61 62 6} 64 

3333 



1111 



65 68 07 63 

3333 



1111 



222 

69 70 71 72 

3333 



1111 

2 222 

33 14 75 76 

3 333 



1111 
2222 

77 73 79 80 

3333 



n+1 m+1 



Punch ed out f r om 2nd Reading St. Punched out from Storage and CR card m+1 

X punched from D.^. 




0000 

12 3 4 

1111 
2222 

12 3 4 

3 3 33 
4444 

12 3 4 

55 5 5 
6666 

12 3 4 

7777 



I 2 3 ft 5 6 7 8 

9 9 9 5j3 y 3 3 3 3 

1 2 3 4'S C 7 E ! 3 10 
I3M 733727 



Fig. 6 Card Diagram of Type 604 Operations 
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of the quantities u and e~( x+y 'n+l, m from the card to general storage 1,2 and fac- 
tor storage 3,4. The transfer goes through the normal side of punch selectors 1,2, 
4, 5, because there is no X-punch in the blank card m+1 in the first reading station. 
During the second calculating cycle, program steps 1-22 are suppressed 
through the normal side of calculator selectors 1,2 (see the calculate control panel 
wiring, Figure 7 ). The operations performed on program steps 23 - 40 are shown in 
the flow diagram. 

A check quantity A ="u - u is the difference between the back solu- 
n J n, m n,m 

tion stored in the counter since the first calculating cycle and the original value of 

u . The two values shall only differ by some units in the last position due to round- 
n,m 

ing errors. As the difference can be either positive or negative, its absolute value 
I A I is entered into the counter by reading a out to general storage 1,2 and subtract- 
ing it into the counter if it is negative. If it is positive, the subtraction is suppressed. 
The absolute value of the difference is tested for zeroes in the second position and up- 
wards. These operations are performed in program steps 23-26. 

The rest of the second calculating cycle is used for computation of a new 
value of e"( x+y 'n+l,m+l according to (6.21) and of n+1 m+1 according to (6.13). Dur- 
ing the second card cycle, e~ Ay • 10 = 99005 was read into the multiplier quotient 
unit from the digit emitter through the normal side of punch selector 3. Calculating 
the new value of e"^ x+y ^n+l, m+1 takes just one multiplication. When this 8-digit num- 
ber is multiplied by u - * to give 0, it has to be split into two parts, called e^ and 
e„. The products 0., and 0„ of e^^ and e 2 ,with u n+1 m+1 appropriately shifted, are ad- 
ded together on program step 39 to make ^ m +l* 

When shifting out the three most significant digits of the product 1 . 10 , 
the validity of (6.18), i.e. , < 1, is checked on program step 36 by a zero test on the 
digits to the left of the decimal point. - ^ is punched from the counter, e~^ x+y 'n+l, m+1 
from general storage 1,2, and u ^ j from general storage 3,4. The punching oc- 
curs during card cycle 1, which belongs to the next computation step. This means that 
the loop in the flow diagram has been closed, which is indicated by the symbol m+1— »m. 

7.4 Thus the calculation proceeds until all the cards in the feed hopper are fin- 

ished. After the 604 run, the old master cards are sorted away for tabulation while 
the new master cards, i.e. , the previous blanks, are interspersed with new blanks. 
The merged deck is preceded by a starting card, which carries an X-punch in column 
70 and the starting value of e~' x+y 'n+2,0 in columns 34-41. The role of the starting 
card, called card 0, is shown in the flow diagram. On the first reading cycle, the X- 
punch picks up calculator selector 1 and punch suppression. In the following calcula- 
ting cycle, program step 1 is the only one which is not suppressed, due to the wiring 
of calculator selectors 1 and 2, as shown in the calculator control panel wiring dia- 
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gram, Figure 7. Program step 1 resets a 2 n and u n+2 to zero in factor storage 
1,2 and general storage 3, 4. Then everything is ready for the first computing step 
of the new integration column. 

The starting card is not punched because of the punch suppression, but it 
is treated as a former blank when reaching the second reading station. As nothing 
has been key-punched in the columns corresponding to u . , this quantity is correct- 
ly set as zero when read into general storage 1,2. The starting value of the exponen- 
tial function, e ~^ x+y 'n+2,0 is read into factor storage 3,4. 

7.5 From the card diagram, Figure 6, it is evident that several of the fields are 

used only when the card is run through the machine for the first time, not the second 
time when it is a master card. This means that the starting deck, representing the 
values on the y-axis, does not need the quantities enclosed in brackets on the master 
cards. The only information necessary to key-punch is n = 0, m and ^/(y) = u n = e m Ay -1. 
The starting deck for the first column is to be distinguished from the starting cards of 
which one has to precede every merged deck of masters and blanks. 



7.6 The programmed checks insure a completely error-free computation and 

punching of u under the condition that the original values of u n = i//(y) are cor- 

' '-PYx+v^ 

rect, which can be checked once and for all. The computation of e KArry ' is easily 

checked on the tabulator by totalling 

1 m=l 

on the same run on which the results are listed. For instance, a mistake in choosing 

the starting card might not influence the computation of u . The error in e~^ x+y 'n,m 

n, m ' 

immediately shows up when checking the total. 

For M-300, S 1 e IlAx has shown to deviate from the exact value l-e' 3 /e°' 01 -l 

by an amount which is almost entirely due to the approximation of e with 5 figures 

in the 604 setup. M M 

The totals of S n , m and jp t n+1, m are listed in the same tabula- 
m=l ' m=l ' 

tion run. Because every such total is listed twice, this will render a further, although 
not necessary, check of the punching of u . The final calculation of <j> according to 
(6. 13) is unchecked due to the shortage of program steps. On a 60-program step mach- 
ine a check could be easily made part of the program. 

8.1 The decimal point convention chosen for u and a sets the limits u < 1000, 

a < 100000. This would do for integration up to (3. 7, 3. 7) in the xy -plane. 

The test integration was made only to y = 3 and x = 0. 43. A check on the ac- 
curacy of the results was obtained through the expansion formula (6. 16). Some of the 
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values are shown in the table below. 



y= 3.00 


Calc. value 


Exact value 


Difference 

(u -z) . io 5 
v n,m ' iU 


X 


u 
n, m 


z(x,y) 


0.01 


19.24698 


19.24697 


1 


0.06 


20.07180 


20.07178 


2 


0.11 


20. 92666 


20.92664 


2 


0.16 


21.81250 


21.81246 


4 


0.21 


22.73025 


22.73021 


4 


0.26 


23.68090 


23.68085 


5 


0.31 


24.66545 


24.66539 


6 


0.36 


25.68493 


25.68486 


7 


0.41 


26.74039 


26.74030 


9 



Table 1 

The accuracy of the calculated exponential function is shown in Table 2 for 
the x-column n=25. The starting value e" n Ax = e' 0, 25 = 0. 7788008. 



x= 0.25 


Calculated 
value 


Exact value 


Difference 


Relative 
Difference 


m 


e -n AX^-niAy 


e -(x+y) 


A.10 7 


A r .10 6 





0.7788008 


0.7788008 






50 


0.4723706 


0.4723666 


40 


8 


100 


0.2865098 


0.2865048 


50 


18 


150 


0.1737787 


0.1737739 


48 


28 


200 


0.1054027 


0.1053992 


35 


33 


250 


0.0639304 


0.0639280 


24 


38 


300 


0.0387756 


0.0387740 


16 


41 



Table 2 



8.2 Although in the test setup the initial function x(x) along the x-axis has been 

chosen as zero, only a minor change is necessary to run the same setup with any re- 
gular function x( x ). 

In the starting card for a new column need only be key-punched the values 



and 



n+1,0 



X( (n+1) A x) 



n,0 



_2_ X '( (n+1/2) A x ) 

Ay 



These values are read into general storage 3,4 and factor storage 1,2 on the starting 
card read cycle. The starting card calculating cycle is skipped. 

Consider the general case of equation (1.11). The only difference from the 
described setup is that an exponential function e -bKx/(ab-c)-ay/K has to be calcula _ 
ted with an equivalent of (6. 21) by multiplying with the constant e~ Ay / K in each step. 
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This means only that other digits are read from the digit emitter into the multiplier 

unit during the second card cycle and that some readjustments of the number of shift 

steps in the following operations might be necessary. The decimal point convention 

for u __ and a „ may also have to be changed, 
n, m n, m J ° 

The calculation of the initial values of x(x) and i// (y) from the correspond- 
ing values of is made once and for all before key-punching the starting decks and 
does not concern the machine setup. 

If it is accepted that the wanted quantity is given as a function of x and y 
instead of its original independent variables | and 17 , it has been demonstrated how 
the canonical homogeneous case of the linear hyperbolic equation with constant coef- 
ficients and characteristic initial values can be set up for numerical integration on the 
604 Electronic Calculating Punch. 
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AN APPROACH TO STORED PROGRAMMING 
ON THE 607 ELECTRONIC CALCULATOR 

Eric V. Hankam 

Watson Scientific Computing Laboratory 
International Business Machines Corporation 



1. Introduction 

The general purpose computing system described in the following paragraphs 
has been designed to make maximum use of the high-speed program unit and elec- 
tronic storage of the 607. In this system as many as six three-address instruc- 
tions can be read from any given card into the 607 storage. The 607 then inter- 
prets and executes these instructions sequentially before the next card is fed. 
Arithmetical, logical, and shift operations can be performed on 8-digit numbers 
(fixed decimal) with sign. It is possible to execute alternative instructions, skip 
instructions, and modify instructions, depending on the result of a logical opera- 
tion. For convenient use of this system, a simple coding scheme has been devel- 
oped. Sections 2 to 5 contain sufficient information to enable a coder to apply the 
system to a problem; sections 6 to 10 are directed toward the reader who is inter- 
ested in the construction of the control panels. 

2. Data 



The data used in a calculation are 8-digit numbers with sign. Up to ten 8-digit 
numbers may be stored in the machine; the ten storage locations are identified by 
the addresses 0, 1,2, 3,4, 5, 6, 7, 8, 9. A number read out of storage is retained in 
storage, and it may be read out again if required. A number read into storage 
will replace the number previously occupying that storage location; i. e. , a storage 
location is cleared automatically prior to entry of a new number. 

Input: Data may be read from cards into storage locations 0, 1,2, 3. 
Output: Results may be punched into cards from storage locations 8, 9. 

Provision has been made to read additional data into storage locations 4, 5, 6 
from a special data card (see sec. 5b). To permit independent entry of data, card 
reading is under the control of X-punches. Similarly, punching is governed by in- 
dividual X-punches in the card fields associated with storage locations 8 and 9 (see 
sec. 5a). It should be noted that a number that is punched from storage will also 
be retained in storage. 
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3. Instructions 

An instruction (A0B=C) is a 5-digit number consisting of three single-digit ad- 
dresses (A, B, C) and a 2-digit operation code 0. The meaning of the A-, B-, and 
C-addresses can be best explained by an illustration. For example, the instruction 
1008 (here, is an arithmetical operation) will cause the machine to search for the 
contents of storage locations 1 and 0, to perform an operation on them, and to read 
the result into storage location 8. In another example, where © represents multi- 
plication, the instruction 3033 will replace the number that is in storage location 3 
by its square. 

Five instructions may be read from a single card. These instructions, which 
are punched into five consecutive card fields, are loaded into 607 storage and execu- 
ted sequentially, i. e. , from left to right as they appear on the card. The number of 
instructions that will be executed from a given card depends on the 2-digit number, 
n, punched in that card. For n=00 (or blank), no instruction is executed. The first 
instruction is always destroyed in the process of execution, but the remaining in- 
structions are circulated in the machine so that, for n >5, instructions 2, 3,4, 5 are 
repeated in a "loop. " A sixth instruction may be included in that loop as explained 
in the next paragraph. After the last instruction has been executed, the next card 
feeds. 

Input and output data may be punched into the card that contains the instructions 
(see card form in sec. 5a). However, an instruction may apply to data already in 
storage as well as to the data read from this card. It is possible to punch a sixth in- 
struction into the card field normally reserved for input data associated with storage 
location 0. 

4. Operations 

a. List of operations and codes 

{00 allow for carry 
Subtract A-B (01 do not allow for carry 

Add A+R i 10 aUow for carr y 

aqq a+b j n dQ not allQw {or carry 

MultiDlv A-B i 20 allow for 16-digit product 

Multiply A a j 21 aUow for 15 _ digit product 

Divide A/B I 30 allowfor IBI ^ IA| < 10 |B| 

Divide A/13 | 31 allowfor | B j>|A| 

Shift T eft A i 40 1 P° sition 

Shift Left A f 41 2 positions 



Balance Test A 



(50 alternative instructions 



(51 modify instructions 
Zero Test A L.. .7. <• -in t skip instructions 
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b. Detailed description of operations 

Subtract A-B — »C, code 00 

The number in the storage location specified by the B-address is 
subtracted from the number in the storage location specified by the 
A-address; allowance is made for a carry, the right-most digit is 
dropped, and the result is read into the storage location specified 
by the C -address. 

xxxxxxxx A 
xxxxxxxx B 




Subtract A-B — >C, code 01 

The number in the storage location specified by the B-address is 
subtracted from the number in the storage location specified by 
the A-address; no allowance is made for a carry, and the result 
is read into the storage location specified by the C-address. 

xxxxxxxx A 
xxxxxxxx B 



xxxxxxxx C 



Add A+B — »C, code 10 

Similar to operation, 00. 



xxxxxxxx A 
xxxxxxxx B 
xxxxxxxxk 



Add A+B — >C, code 11 

Similar to operation, 01. 

xxxxxxxx A 
xxxxxxxx B 
xxxxxxxx C 

Multiply A- B— >C, code 20 

Multiply A by B, allow for a 16-digit product, round and drop 8 
digits, and store product in C. 

(xxxxxxxx) (xxxxxxxx) = xxxxxxxx xxxxxxxx 



B C 5 



^round 
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Multiply A-B — *C, code 21 

Multiply A by B, allow for a 15-digit product, round and drop 7 
digits, and store product in C. 

(xxxxxxxx) (xxxxxxxx) = Oxxxxxxxx 

A B C -. 

round 

Divide A/B — *C, code 30 

Divide A by B, allow for 9 quotient digits, round and drop rightmost 

digit, and store quotient in C. 

If |B|i |A| < 10 |b|, the quotient will have 8 significant digits. 

If |B|> |A|, the quotient will have <8 significant digits. 

This operation will cause a quotient overflow and result in a zero 

quotient for all cases | Af =10 |B| 

A xxxxxxxx- 10 8 = XXXXXXX2 ^ 

B XXXXXXXX ^ " ^A 

*roimd 

Divide A/B— *C, code 31 

Divide A by B, allow for 9 quotient digits, round and drop rightmost 
digit, and store quotient in C. 

If |A| < B g 10 |A|, the quotient will have 8 significant digits. 
If |B| > 10 |A|, the quotient will have <8 significant digits. 
This operation will cause a quotient overflow and result in a zero 

quotient for all cases | A| £ |B| . 

9 
A xxxxxxxx ♦ 10 = Oxxxxxxxxk 

B xxxxxxxx a, r 

u Ground 

Shift Left A— >C, code 40 

A is shifted one place to the left and the result is stored in C. 
Using a C-address equal to the A-address has the effect of return- 
ing the shifted number to its original storage location. The B-ad- 
dress has no function in this operation and can be left blank. 

Shift Left A— >C, code 41 

A is shifted two places to the left, and the result is stored in C. 
As in operation 40, the B-address has no function in this operation 
and can be left blank. 
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Balance Test A, code 50 

The algebraic sign of A is tested. If A is negative, instructions 1 and 2 
will be ignored in subsequent cards until a card that has a Y-punch in 
column 5 is read; instructions 1 and 2 in the Y^-card will also be ignored, 
but, thereafter, they will be executed in the normal manner. For example, 
when n=2, this "branching" operation will cause instructions 1 and 2 to be 
executed if A is positive, or instructions 3 and 4 if A is negative. 
As a by-product of the design of the control panel, B will transfer into C. 
In general, both the B-address and C-address will be left blank. 

Balance Test A, code 51 

The algebraic sign of A is tested. If A is negative, the A-, B-, and C-ad- 
dresses of instructions in all subsequent cards will be modified by a. , , 
and y, respectively, until a card with a Y-punch in column 6 is read. In- 
structions in the Y g -cards will also be modified, but, thereafter, they will 
be executed in the normal manner. The single digit numbers, a , /3, and y, 
are punched in the card and are applied additively to the addresses of an 
instruction to transform them to A+a , B+/3 , and C+ y . The section on 
the 607 program shows that the addition is performed as follows: 





C 


I 


> 


B 


A 




7 









a 


X 


X 


x 


x 


X 


X 



A carry resulting from (A+a) will increase (B+/3 ); likewise © may be in- 
creased by 1. If C + y > 9, it will be treated as 9. By means of a single 
X-punch, minus signs may be attached to a, , and y . Though 
a , |3 , y may change from card to card, all instructions in a single 
card will be modified by the same at , fi , y . 

As a by-product of the design of the control panel, B will transfer into C. 
In general, both the B-address and C-address will be left blank. 

Zero Test A, code 60 

A is tested for zero. If A is zero, all subsequent instructions in the same 
and following cards are skipped until a card with a Y-punch in column 7 is 
read. Instructions in the Y„-card will also be skipped, but, thereafter, 
they will be executed in the normal manner. 

As a by-product of the design of the control panel, B will transfer into C. 
In general, both the B-address and the C-address will be left blank. 
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Zero Test A, code 61 

The only difference between this operation and operation 60 is that all posi- 
tions except the rightmost position of A are tested for zero, i.e. , opera- 
tion 61 tests whether |A| < 10. 

c. Further remarks on operations. 

A 7-, 8- or 9-punch in the tens position of the operation code will be treat- 
ed as "No Operation. " In this operation the A-address has no function, but B 
may be transferred into C. A number may also be transferred from one stor- 
age location to another by the addition of zero to it, but the former method has 
the advantage of not requiring a storage location for the number zero. 

The sign of a number can be changed by the subtraction of the number from 
zero. 

A number may be shifted to the right, one place at a time, by the addition 
of zero to it by means of operation 10. Since the low-order positions of a prod- 
uct are dropped, we may shift a number k places to the right by multiplying it 

8-k 
by 10 (where k = 1, 2, . . . , 8) with operation 20. In the latter case, rounding 

will take place if the right shift causes significant digits to be dropped. 

Since the 607 reads blank columns as zeros, an unpunched position in an 
instruction will always be treated as a zero. In particular, the machine will 
consider a completely blank instruction as 00 0, and if it executes such an 
instruction (subject to n), it will, in effect, clear storage location 0. As a mat- 
ter of fact, any storage location can be cleared by specification of its address 
as the C-address in an otherwise blank instruction. 

A study of the section describing the 607 program will show that the mach- 
ine simply tests whether the units position of the operation code is zero or not; 
therefore, the 1 in the units position is really arbitrary, and any other non-zero 
digit can be substituted. 

Card Form 

a. General card form for instructions and data 



a 
o 

+j 

a 
o 
.^ 

•.-1 
-M 

c 


n 


YYY 


instr 

#1 


instr 

#2 


instr 


instr 

M 


instr 

#5 


XX 

data 
into 
storage 




X 

data 
into 
storage 

1 


X 

data 
into 
storage 

2 


X 

data 
into 
storage 

3 


X 

punch 
from 
storage 

8 


X 
punch 

from 
storage 

9 


T3 

1-2 


3-4 


5-7 


8-12 


13-17 


18-22 


23-2*5 


28-31 


f 33-40 


41-48 


49-56 


57-64 


65-72 


73-80 
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Col. 1,2 : identification 

Col. 3, 4 : n (number of instructions to be executed) 

Col. 5-7 : a , /3 , y , (address modifiers). For a complete explanation 
of a , p , y , and the significance of the Y-punches see the 
description of the Balance Test and Zero Test operation (sec. 4b). 
An X- punch in col. 7 will make a ,$, y negative. 

Col. 8-32 : five 5-posit ion instruction fields. The leftmost position of each 
field contains the A-address; the 2nd and 3rd positions contain 
the operation code 0; the 4th position contains the B-address; 
and the last position contains the C-address. 

Col. 33-64 : four 8-position input-data fields are used to enter data into 
storage locations 0, 1,2,3, respectively. An X-punch in the 
leftmost position of each field will cause the corresponding 
storage location to read in. An X-punch in the rightmost posi- 
tion of the field identifies negative numbers. 
An X-punch in col. 35 makes it possible to use col. 36-40 for 
a sixth instruction. An X-punch in col. 33 is also required to 
read the instruction into storage location 0. 

Col. 65-80 : two 8-position output-data fields are used to punch results from 
storage locations 8,9, respectively. An X-punch in the leftmost 
position of each field will cause the corresponding storage loca- 
tion to read out. An X-punch in the rightmost position of the 
field identifies negative numbers. 

b. Special data card form ( Y in col. 1) 



/Y 


x 


X 


X 


/ 




data 
into 
storage 
4 


data 
into 
storage 
5 


data 
into 

storage 
6 




1 


9-16 


17-24 


25-32 


/ 



A Y-punch in col. 1 makes it possible to read data from three adjacent 
8-position fields (col. 9-32) into storage locations 4, 5,6, respectively. An 
X-punch in the leftmost position of each field will cause the corresponding 
storage location to read in. An X-punch in the rightmost position of the field 
identifies negative numbers. 

Instructions cannot be read from this card, but instructions already in 
storage will not be destroyed, and can be executed if n * is punched into 
col. 3,4. 



75 



c. Further remarks on card forms. 

Card forms are quite arbitrary, and we can alter them readily by 
making corresponding changes on the 529 control panel. Suggestions for some 
modifications appear in the Appendix. 

607 Storage 

a. Data Storage 

The following 607 storage units have been used for storage loca- 
tions 0,1, . . . , 9: 

0~ Gl,2 

1 ~ G3,4 

2 ~ G9,10 

3 ~ Gil, 12 
4~ G13,14 

5 ~ G15,16 

6 ~ G17,18 

7 ~ G19,20 

8 ~ G21,22 

9 ~ G23,24 

The 607 at the Watson Scientific Computing Laboratory is a model 
with 140 program steps and, therefore, does not have G5,6,7, 8. 

b. Working Storage 

Fl,2 stores A 

F3, 4 stores B before an operation and C after an operation 

MQ 

EC 

c. Miscellaneous Storage 

F5 stores n Two card columns have been reserved for n. 

If a third card column were included on the 529 
control panel, the 607 program could handle any 

F7 stores a@y n g 999. The 529 control panel has been wired 
so that a and y are interchanged as they are 
read from the card into F7, i. e. , 
on card 



a <3 y 



in F7 
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F6 


instruction 1 


F8 


2 


F10 


3 


F12 


4 


F9,ll 


5 



d. Instruction storage 

The 529 control panel has been wired so 
that every instruction is read from the 
card into storage as follows: 

A © B C on card 

C B A in storage 

Gl,2 may also be used for instruction 6. 

Explanation of the 607 program 

In the following, the ith program step will be denoted by P. . 

a. Pj to Pg and P. ~: n is reduced by 1 and tested. If n g 0, program re- 
peat is impulsed; if n <0, P 1Q will be active and terminate the program. 

b. P 4 to P 9 and P n : P 4 to P ? transform y a to y 00 /3 a, and on P g the 
instruction in F6 is added to it. As mentioned in section 6d, all instructions 
in storage are of the form C B A, and it will be shown later (sec. 7h) that 
every instruction is put into F6 prior to its execution. P g and P u will store 
the C0B part of the instruction in F6, and leave the isolated A-address in EC 
(electronic counter). 

c. P 12 to Pg Q : Interpretation of the A-address: By successive subtractions 
of l's and balance tests, the machine interprets the A-address, reads out the 
appropriate general storage unit, and puts A into Fl, 2. Once A has been found, 
subsequent program steps (in the group P 12 to P 3Q ) are suppressed with the 
aid of GS # 1 (Group Suppression # 1). 

d. P 31 to P 52 : Interpretation of the B-address: P gi to P, 3 put the B-ad- 
dress into EC and leave C0, the remainder of the instruction, in F6. Then 
(P 34 to P 52 ) by successive subtractions of l's and balance tests, the machine 
interprets the B-address, reads out the appropriate general storage unit, and 
puts B into F3, 4. Once B has been found, subsequent program steps (in the 
group P 34 to P 52 ) are suppressed with the aid of GS# 1. 

e. Pg 3 to Pgg : The machine interprets the units position of the operation 
code. The units position is zero tested (Prg), and the result of the test is ap- 
plied to later steps. 

f. Pg^ to Pj 14 : The machine interprets the tens position of the operation 
code, performs the appropriate operation, and, except for logical operations, 
puts the result (C) into F3, 4. 
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a . Pc 7 to P 7Q : The tens position of the operation code is placed 
into EC (P 57 to P 59 ). If the tens position is 0, the operation is 
subtraction, and the sign of B is changed so that the addition rou- 
tine can be used later for both addition and subtraction (Pg Q to Pg 3 ). 
If it is 1 or 0, the machine adds Fl,2 to F3,4 and puts the result 
(subject to the zero test on P,-g) into F3,4 (P g4 to P 70 )- 
/3. P 71 to P g4 : If the tens position of the operation code is 2, the 
machine multiplies Fl,2 by F3,4. The product is rounded (subject 
to the zero test on Pgg) and read into F3,4. 

y. P 8 c to Pg 8 : If the tens position of the operation code is 3, the 
machine divides Fl, 2 by F3, 4. The positioning of the dividend is 
subject to the zero test on P 5g . The quotient is rounded and read 
into F3,4. 

6. Pg 9 to P 1Q2 : If the tens position of the operation code is 4, the 
number in Fl,2 is transferred into F3,4 and at the same time shifted 
one or two places to the left (subject to the zero test on Peg). 
e# P 10 „ to P 1Q7 : If the tens position of the operation code is 5, the 
sign of the number in Fl, 2 is tested. If negative and the units posi- 
tion of the operation code was zero, N.B. selector 1 is picked up. 
If negative and the units position of the operation code was not zero, 
N.B. selector 2 is picked up.* 

£• P 1Q8 to P 114 : K the tens P° sition °* tne operation code is 6, the 
machine tests whether the number in Fl,2 is zero. Subject to the 
zero test on P 56 , either the entire number in Fl,2 or the entire num- 
ber except the rightmost position is tested. If zero, N.B. selector 
3 is picked up.* Furthermore, F5 is cleared (P 114 ) which will cause 
any remaining instructions to be skipped. 

g. P 115 to P 134 : Interpretation of the C-address: By successive subtract- 
ions of l's and balance tests, the machine interprets the C-address and reads 
F3,4 (C) into the appropriate general storage unit. Once this storage unit has 
been determined, subsequent program steps (in the group P 115 to P 134 ) are 
suppressed with the aid of GS # 1. 



* The use of N. B. selectors is explained in the section on the 529 control 
panel (see sec. 9b). 
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h " P 135 to P 140 : The instruction s are "circulated" within the machine by 
advancing each instruction into the storage unit occupied by its predecessor. 
The instruction to be executed on the next program sweep is put into F6; this 
instruction is also moved into F9, 11 (if calculator selector 1 is normal) so that 
it can be reused (modulo 4). Calculator selector 1 can be picked up on the 529 
control panel by an X-punch in col. 35 (see sec. 9b), and if it is transferred, 
the instruction in F6 is moved into Gl,2 so that it can be reused (modulo 5). 

Another comment concerning this group of steps will be found in the 
Appendix. 



8. The 607 Program * 



1: 

2: 

3: 

4: 

5: 

6: 

7: 

8: 

9: 

10: 

11: 

12: 

13: 

14: 

15: 

16: 

17: 

18: 

19: 

20: 

21: 

22: 

23: 

24: 

25: 

26: 

27: 

28: 

29: 

30: 

31: 

32: 

33: 

34: 

35: 

36: 

37: 

38: 

39: 

40: 



F5RO; 
Emit 1; 
EC RR; 



F7 

EC 

MQ 

MQ 

F6 

EC 



RO; 

RO; 

RO; 

RO; 

RO; 

RO; 
Prog. End; 
F6 RO; 
Emit 1; 
Gl,2RO; 
Emit 1; 
G3, 4 RO; 
Emit 1; 
G9, 10 RO; 
Emit 1; 
Gll,12RO 
Emit 1; 
G13, 14 RO 
Emit 1; 
G15,16 RO 
Emit 1; 
G17,18RO 
Emit 1; 
G19,20RO 
Emit 1 ; 
G21,22 RO 
G23,24RO 
F6 RO; 
EC RO; 
F6 RO; 
Emit 1 
Gl,2RO; 
Emit 1 
G3,4RO; 
Emit 1 
G9,10RO; 
Emit 1 



EC RI+ 

EC RI- ; 

F5 RI; 

EC RI+ ; 

MQ RI; 

EC RI- ; 

EC RI+ ; 

EC RI+ ; 

F6 RI; 

EC RI- ; 
EC RI- ; 
Fl,2 RI; 
EC RI- ; 
Fl,2 RI; 
EC RI- ; 
Fl,2 RI; 
EC RI- ; 
Fl,2 RI; 
EC RI- ; 
Fl,2 RI; 
EC RI- ; 
Fl,2 RI; 
EC RI- ; 
PI, 2 RI; 
EC RI- ; 
Fl,2 RI; 
EC RI- ; 
Fl,2 RI; 
Fl,2 RI; 
EC RI+ ; 
F6 RI; 
EC RI- ; 
EC RI- ; 
F3,4RI; 
EC RI- ; 
F3,4RI; 
EC RI- ; 
F3,4RI; 
EC RI-; 



B.T. 
Prog. Rep. 



r'e 



Ro 3rd 
Ri 3rd 
Ri 5th 

Ro 2nd 

Ri 2nd 
B.T. ; 
EC RE 
B.T. ; 
EC RE 
B.T 
EC 
B.T.; 
EC RE 
B.T. ; 
EC RE 
B.T. ; 
EC RE 
B.T.; 
EC RE 
B.T. ; 
EC RE 
B. T. ; 
EC RE 
EC RE 

Ro 2nd 
Ri 2nd 
B.T. ; 
EC RE 
B. T.; 
EC RE 
B.T.; 
EC RE 
B.T. : 



;Sup- 
Sup- 
Sup- 
Sup- 
Sup- 
Sup- 
Sup- 
Sup+ 

GS#1 

GS#1 

Sup- 

GS#1 

Sup- 

GS#1 

Sup- 

GS#1 

Sup- 

GS#1 

Sup- 

GS#1 

Sup- 

GS#1 

Sup- 

GS#1 

Sup- 

GS#1 



d.o. 
p.u. 

p.u. 

p.u. 

p.u. 

p.u. 

p.u. 

p.u. 

p.u. 

p.u. 



GS#1 d.o. 

GS#1 p.u. 

Sup- 

GS#1 

Sup- 

GS#1 

Sup- 



p.u. 



p.u. 



Sup+ 
Sup+ 
Sup+ 
Sup+ 
Sup+ 
Sup+ 
Sup+ 
Sup+ 
Sup+ 



Sup#l 

Sup#l 

Sup#l 

Sup #1 

Sup# 1 

Sup# 1 

Sup# 1 

Sup# 1 
Sup# 1 



Sup+ 

Sup+; Sup # 1 

Sup+; Sup # 1 



The notation in this and following sections is similar to that employed in Applied 
Science Newsletter No. 3, pp. 41-42. 
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41: 

42: 

43: 

44: 

45: 

46: 

47: 

48: 

49: 

50: 

51: 

52: 

53: 

54: 

55: 

56: 

57: 

58: 

59: 

60: 

61: 

62: 

63: 

64: 

65: 

66: 

67: 

68: 

69: 

70: 

71: 

72: 

73: 

74: 

75: 

76: 

77: 

78: 

79: 

80: 

81: 

82: 

83: 

84: 

85: 

86: 

87: 

88: 

89: 

90: 

91: 

92: 

93: 

94: 

95: 

96: 

97: 

98: 



Gil, 12 RO 
Emit 1; 
G13,14RO 
Emit 1; 
G15,16 RO 
Emit 1; 
G17,18RO 
Emit 1; 
G19,20RO 
Emit 1; 
G21,22RO 
G23, 24 RO 
F6RO; 
EC RO; 
F6RO; 
Z.T.; 
F6RO; 
EC RO; 
F6RO; 
Emit 1; 
EC RE; 
F3,4RO; 
EC RR; 
Emit 1; 
EC RE; 
Fl,2RO; 
F3,4RO; 
EC RR; 
EC RR; 
Skip On; 
Emit 1; 
F3,4RO; 
F3,4 RO; 
EC RR; 
Fl,2 RO; 
F3,4RO; 
EC RR; 
F3,4RO; 
Fl, 2 RO; 
1/2 Adj; 
EC RR; 
1/2 Adj; 
EC RR; 
Skip On; 
Emit 1; 
EC RE; 
Fl,2 RO; 
Fl,2 RO; 
F3,4RO; 
EC RR; 
Fl,2RO; 
MQRO; 
F3,4RO; 
MQRO; 
Fl,2RO; 
1/2 Adj; 
EC RR; 
Skip On; 



F3,4RI; 
EC RI- ; 
F3,4RI; 
EC RI- ; 
F3,4RI; 
EC RI- ; 
F3,4RI; 
EC RI- ; 
F3,4RI; 
EC RI- ; 
F3,4RI; 
F3,4RI; 
EC RI+ ; 
F6RI; 
EC RI- ; 
EC RE 
EC RI+ 
F6RI; 
EC RI- ; 
EC RI- ; 

EC RI- 
F3,4RI 
EC RI- 

ECRI+ 
EC RI+ 
F3,4RI 
F3,4RI 

EC RI- ; 

MQRI; 

EC RI+ ; 

F3,4RI; 

MPL+; 

MQRI; 

F3,4RI; 

EC RI+ 

MPL+; 



F3,4RI; 

F3,4RI; 

EC RI- ; 

EC RI+ ; 
EC RI+ ; 
DIV. 
F1,2RI; 
EC RI+; 
Fl,2 RI; 
DIV; 
EC RI+; 
EC RI+ ; 

F3,4RI; 



EC RE 
B.T.; 
EC RE 
B.T.; 
EC RE 
B.T.; 
EC RE 
B.T.; 
EC RE 
B.T.; 
EC RE 
EC RE 

Ro 2nd 
Ri 2nd 



Ro 2nd 
Ri 2nd 
B.T.; 



B.T. ; 

Ro 2nd 

B.T. 
EC RE; 

Ro 6th; 
Ro 6th 



Ri 3rd 
Ro4th 
Ri 2nd; 
Ro 3rd 

B.T. 

Ri 5th 
Ri 6th 



Ri 5th 
EC RE; 
Ri 5th 
Ro 2nd; 



GS # 1 p.u. 

Sup- 

GS # 1 p.u. 

Sup- 

GS # 1 p.u. 

Sup- 

GS # 1 p.u. 

Sup- 

GS # 1 p.u. 

Sup- 

GS # 1 p.u. 



Supf 
Sup+ 
Sup+ 
Sup+ 
Sup+ 
Supf 



Sup #1 

Sup #1 

Sup #1 

Sup # 1 

Sup #1 

Sup# 1 
Sup# 1 



Supf 
Sup+ 
Supf 

Supf 

Supf 

Supf 

Sup+; 

Sup+; 

Supf 

Sup+ 

Sup+ 

Supf 

Sup+ 

Sup+ 

Sup+ 

Supf 

Sup+ 

Sup+; 

Supf; 

Supf; 

Sup+; 

Supf 

Sup+ 

Supf; 

Sup+; 

Sup+ 

Supf 

Sup+ 

Sup+ 

Sup+ 

Supf 

Supf 

Supf 

Supf 

Supf 



SupNZ 
Sup Z 



SupNZ 
SupNZ 
Sup Z 
Sup Z 



SupNZ 
Sup Z 
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99: Emit 1; 
100: Fl,2 RO; 
101: Fl,2 RO; 
102: Skip On; 
103: Emit 1; 
104: EC RE; 
105: Fl,2 RO; 
106: Fl,2 RO; 
107: Skip On; 
108: Emit 1; 
109: Fl,2 RO; 
110: EC RR; 
111: Fl,2 RO; 
112: Z.T. ; 
113: Emit 1; 
114: F5 RI; 
115: F6 RO; 
116: Emit 1; 
117: F3,4RO 
118: Emit 1; 
119: F3,4RO 
120: Emit 1; 
121: F3,4 RO 
122: Emit 1; 
123: F3,4RO; 
124: Emit 1; 
125: F3,4RO; 
126: Emit 1; 
127: F3,4RO; 
128: Emit 1; 
129: F3,4 RO; 
130: Emit 1; 
131: F3,4RO; 
132: Emit 1; 
133: F3,4RO; 
134: F3,4RO; 
135: F8 RO; 
136: F10 RO; 
137: F12 RO; 
138: F9,ll RO; 
139: Gl, 2 RO; 
140: F6 RO; 



EC RI- ; 
F3,4RI; 
F3,4RI; 

EC RI- ; 

EC RI+ ; 
EC RI+ ; 

EC RI- ; 

EC RI+ ; 

Fl,2 RI; 

EC RI+ ; 

EC RE; 

EC RI- ; 

EC RI+ ; 
EC RI- ; 
G1,2RI; 
EC RI- ; 
G3,4RI; 
EC RI- ; 
G9, 10 RI; 
EC RI- ; 
G11,12RI; 
EC RI- ; 
G13,14RI; 
EC RI- ; 
G15,16RI; 
ECRI- 
G17, 18 RI; 
EC RI- ; 
G19,20RI; 
EC RI- ; 
G21,22 RI; 
G23,24RI; 
F6RI 
F8 RI 
F10RI 
F12 RI 
F9, 11 RI; 
F9, 11 RI (Cal-Sel 1 N); Gl', 2 RI (Cal-Sel 1 T) 



B.T. 
Ri 2nd; 
Ri 3rd; 

B.T. 

EC RE; 
EC RE; 

B.T. 

Ro 2nd; 



EC RE 
Skip Off 
B.T.; 
EC RE 
B.T.; 
EC RE 
B.T.; 
EC RE 
B.T.; 
EC RE 
B.T.; 
EC RE 
B.T.; 
EC RE 
B.T.; 
EC RE 
B.T.; 
EC RE 
B.T.; 
EC RE 
EC RE 



Sup+; 
Sup*-; 
Sup+ 

Sup+ 

Supf; 

Sup+; 

Sup+ 

EC RE 

Sup+ 

Sup+; 

Sup+; 

Sup+ 

Sup+; 

Sup+; 

GS#1 d.o. 

GS# 1 p.u. 

Sup- 

GS#1 

Sup 

GS#1 p.u. 

Sup- 

GS#1 

Sup- 

GS#1 

Sup 

GS#1 

Sup 

GS#1 

Sup- 

GS#1 

Sup- 

GS#1 



p.u. 



p.u. 



p.u. 



p.u. 



p.u. 



p.u. 



p.u. 



Sup NZ; 
Sup Z; 



Sup NZ; 
Sup Z; 



Sup Z 
Sup Z 

Sup NZ; 
Sup NZ 



Sup+ 

Sup+; 

Sup+; 

Sup+; 

Sup+; 

Sup+; 

Sup+; 

Sup+; 

Sup+; 



Sup w/o Test (Cal-Sel 1 N) 



EC RE 
EC RE 



N.B.#1 p.u. 
N.B.#2 p.u. 



N.B.#3 p.u. 



Sup# 1 

Sup# 1 

Sup# 1 

Sup# 1 

Sup# 1 

Sup# 1 

Sup# 1 

Sup# 1 
Sup# 1 



Storage Assignments 

"8-6": Gl, 3, 9,11, 13,15,17,19,21, 23, Fl,3. 

"6-4": F9 
NOTE: Program expanders are used whenever more than four program 
exits are required. 

The 529 Control Panel 



This section outlines the main features of the 529 control panel. 
A review of sections 5 and 6 will show the card fields that are wired to var- 
ious 607 storage entries and exits. A summary is listed below. 
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a. Entries and Exits. 

n >F5* ; C.C. > F5 RI 

a p y >F7* ; C.C. »F7RI 

Instruction 1 — > F6* 

2 — »F8* 

3 — >. F10* 

4 _» F12* 

5 _^F9,11* 



>> 



Instructions are wired through the normal side of a group of selectors 
picked up by Y<, **. The F6, 8, 9, 10, 11, 12 RI hubs are impulsed by C. C. 
also through the normal side of these selectors. Through the transferred 
side of the selectors, three 8-position input-data card fields are wired to 
the entries of G13, 14; G15, 16; and G17, 18, respectively. X-punches in the 
leftmost position of each field are wired to the RI hubs of the corresponding 
storage units (for further details see sec. 5b): 



X 33 >G1,2 RI 

X 41 *G3,4 RI 

X 49 >G9,10 RI 

X 57 > Gil, 12 RI 

G21, 22 ► punching hubs; C. C. — » G21, 22 RO* 

G23, 24 — ¥ punching hubs; C. C. f G23, 24 RO* 



Input-data > Gl,2* 

» >G3,4 

" >>G9,10 

» > Gil, 12 



To clear the counter, wire C.C. — ¥ EC RR 

b. Selectors 

G21,22 RO is wired through the transferred side of a pilot selector 

picked up by Xgg- 
G23, 24 RO is wired through the transferred side of a pilot selector 

picked up by X»g. 
Y* (special data card) is wired to the I-p. u. of a group of punch 

selectors. 
The use of these selectors has already been explained (see sec. 9a). 
X„, is wired to the p.u. hub of calculator selector 1, which permits the 
use of Gl,2 for a sixth instruction. It will be recalled (see sec. 6d) that when 
an instruction is read into storage its A-address and C-address are inter- 
changed; a punch selector, picked up by X 35 , is used to transpose the leftmost 
and rightmost position of the sixth instruction. 



* Through selectors as explained in section 9b. 
** Y. means Y-punch in column i. 
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Latch selectors: As explained in section 7, N. B. selectors 1,2 and 3 
may be picked up as a result of certain negative or zero conditions. Each 
N. B. selector is used to transfer a pilot selector. With the aid of an auxi- 
liary selector, the pilot selector is held in a transferred position until a card 
with a special Y-punch is read. The pilot selector will also be transferred 
for the Y-card, but, thereafter, it will return to its normal position. Punch 
selectors can be coupled to the pilot selector. The diagram below illustrates 
a method for wiring a typical "holding circuit." 

Pil-Sel auxiliary 

X or D-p. u. 9 K Pil-or Pch-Sel 

I -p. u. 



'Y (from card) 
T 
N.B. SeU^-^A N • _^~-» N 

C 

Coupling Exit p.u. '""^hot "X" or Digit 

Pch-Sel 

We shall refer to the group of selectors associated with N. B. selectors 
1,2,3 as latch selectors I, II, III, respectively. Y 5 , Y g , Y„ pick up the aux- 
iliary selectors and initiate the "drop out" of latch selectors I, II, III, re- 
spectively. 

Latch selector I is used to prevent the entry of instructions 1 and 2 into 
F6 and F8. In their place instructions 3 and 4 are entered into F6 and F8. 
(See Appendix for further remarks regarding the selection of instructions. ) 

Latch selector II permits the entry of a/3y into F7, i.e. , a/3y is wired 
through the transferred side of latch selector II. Since the RI hub of F7 is im- 
pulsed by C.C. all the time, F7 will normally be zero. 

Latch selector III prevents the entry of n into F5, i.e. , n is wired 
through the normal side of latch selector III. This is equivalent to reading 
n=0 into F5 whenever latch selector III is transferred, since the RI hub of F5 
is always impulsed by C.C. 
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10. Appendix 

a. Suggestions for modifications. 

Unlike data entry, which is under the control of individual X-punches, 
instructions are read automatically into storage. The card-cycles pulse, 
which is responsible for the automatic entry of all instructions, could be re- 
placed by X-punches to control the entry of each instruction individually. The 
user would have to decide whether the flexibility gained thereby would com- 
pensate for the inconvenience of punching additional X's into the card. 

In the present set-up, only two results may be punched out from G21, 
22 and G23, 24. This number could be increased by the proper selection of 
data input and output. 

The selection of instructions by means of latch selector I need not be 
done entirely on the 529 control panel. Several punch selectors may be elim- 
inated by the association of a single calculator selector, say Cal-Sel 2, with 
latch selector I. Only the selection between instruction 1 (through the normal 
side of latch selector I) and instruction 3 (through the transferred side of 
latch selector I) has to be done on the 529. The selection of the remaining in- 
structions may be accomplished by the following minor modification in the 607 
program. 

P 13? : F12 RO; F10 RI (Cal-Sel 2 N); F6 RI (Cal-Sel 2 T) 

For maximum speed, Sup w/o Test (Cal-Sel 2 T) should be added to Pi 35 
and P 136 . 

b. Suggestions for additional operations and instructions. 

With the aid of calculator selectors, additional operations could be in- 
corporated. Selected program exits can be used for multiple purposes, and an 
effort should be made to design the 607 control panel so that a group of steps 
can be used for more than one operation even without selection. However, since 
calculator selectors stay picked up for the entire calculation, the coder would 
be required to exclude on any given card any instructions that involve both the 
normal and transferred positions of a calculator selector. If selectors which 
could be picked up and dropped out on different program sweeps, such as the 
program repeat selectors on the 605, were available, program exits could be 
selected so as to produce more than 140 distinct program steps. 

The number of instructions that may be read from a single card depends 
primarily on the storage capacity of the machine. Each instruction requires at 
least one 5-position storage unit. In the preceding sections it has been shown 
how a storage unit can be used alternatively for data and instructions. Conse- 
quently, more instructions may be included at the expense of data storage. The 
ease with which additional instructions may be incorporated in the 607 program 
is evident from a study of program steps 135-140. Only a single additional pro- 
gram step is required for each additional instruction. 
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AN ITERATIVE TECHNIQUE FOR MULTIPLE CORRELATION ANALYSIS 

Martin H. Greenberger * 
Joe H. Ward, Jr. 

Air Force Personnel and Training Research Center 
Lackland Air Force Base 
San Antonio, Texas 



I. Introduction 

An earlier iterative method by the same authors appeared in IBM Technical 
Newsletter No. 6 and may be referred to as background for this paper. The same no- 
tation has been maintained as much as possible. One exception is that V. has been re- 
placed by r. Q to represent the ith criterion correlation coefficient. 

Although the approach here is somewhat different, the problem is identical 
with the one treated by the earlier method. It is desired to "fit" (in the least squares 
sense) a hyper-plane in (n+1) space to the set of N points 

(Z 10' Z ll» Z 12> '••' Z ln> 
( Z 20' Z 21' Z 22' *••' Z 2n) 

• • • • • 

(Z N0' Z N1' Z N2' "•' Z Nn) 

in order to permit prediction of the Oth dimension coordinate from the coordinates of 
the other n dimensions. That is, given any arbitrary point (Z ,, Z «,..., Z ), 
a least squares set of weights Bj, B 2 , . . . , B n is desired which will predict Z Q by 
the equation 

Z n = Bj Z - + B Z n + . . . + B Z 

mO 1 ml 2 m2 n mn 

In practice, N will be much larger than n. Hence by replacing the original N points 
with Bj, Bg, • . . , B n , considerable data reduction is obtained. 

Both iterative methods begin with a set of weights X., , X„, . . . , X which 
fail to yield the desired fit. By successive corrections on one X at a time the fit is 
approached. If 



* Now with the Harvard Computation Laboratory. 

Note: The opinions or conclusions contained in this report are those of the authors. 
They are not to be construed as reflecting the views or indorsement of the 
Department of the Air Force. 
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L l = X l Z ll + X 2 Z 12 + ' * ' + X n Z ln 



L 2 = X l Z 21 + X 2 Z 22 + * ' ' + X n Z 2n 



L N =X 1 Z N1 + X 2 Z N2 + '-* + X n Z Nn' 

the desired fit can be achieved by reducing the differences between the Lj and the 
Z . with the appropriate succession of corrections. 

The earlier method converged on the solution by choosing at each stage a 
correction on some selected X which minimized 

N 2 

. S (Z 10 - L.) 2 

with all other X held constant. In contrast, the present method chooses that correc- 
tion for the selected X which maximizes the correlation between (Lp L 2 , . . . , L N ) 
and (Z 1QJ Z 2Q , . . . ,. Z NQ ), all other X fixed. In both cases, the derived least squares 
weights are labeled B,, B„, . „ . , B . 

The former iterative method afforded the utmost in simplicity computation- 
ally, and was therefore ideal for 602 -A application. The present method involves 
more sophisticated computation but offers faster convergence. In theory any number 
of variables may be treated by this technique. The present program handles up to 75 
variables. 

The iterative technique to be developed was designed for the 607 Electronic 
Calculator. However, with some revision it may be adapted to the 650. The present 
approach eliminates the necessity for panel wiring changes during the operation. From 
iteration to iteration there is completely automatic field selection. It has been in suc- 
cessful operation since June 1954, and has provided precious time savings during this 
period. Section VH discusses one of the representative problem solutions. 

The procedure is an offspring of the Kelley-Salisbury method* with some 
important modifications. For example, it is not satisfactory to correct at each stage 
on the equation with the largest numerical error, as the Kelley-Salisbury method sug- 
gests. It is wholly within reason for a variable to produce the greatest error without 
admitting of an appreciable correction by the criterion used. One further modification 
is the punching of e. (defined below) instead of a r The former has the advantage of 
being restricted to a narrow range of magnitude. 



* Kelley, T.L. and Salisbury, F.S. , "An Iteration Method for Determining Multiple 
Correlation Constants, " Journal of the American Statistical Association, 21:282-292, 



1926. 
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It would be well at this point to introduce some additional notation which 
will be used throughout the ensuing development. The following relations are satis- 
fied by the least squares weights: 

B l + r 12 B 2 + • • • + r ln B n = r 10 
r 12 B l + B 2 + . . . + r 2n B n = r 2Q 



r ln B l + r 2n B 2 + -'- + B n = r n0 



where r.. = the correlation between Z. and Z. 
i] i J 

r ij -'jiandPjj-1. 



For the intermediate weights X* , X , . . . , X„ the following a. are defined: 

x Ci n i 



X l + r 12 X 2 + • • • + r ln X n " a l 
r 12 X l + X 2 + . . . + r 2n X n = a 2 



r. X« + r „ X„ + . . . + X„ = a„ 

In 1 2n 2 n n 



Unless Xh , X 2 , . . . , X is identical with B, , B„, , B the a. will not in general 

equal the r.„ . Letting a. - r.,. = e., and calling this difference the ith error, the 
condition for solution is that the e. approximate zero for all i = 1 to n. 

II. Mathematical Statement 

A. We are given the n + 1 variables z n , z., , z„, . . . , z , where Z. = (Z... , Z. , . . . , Z. ) 
is a set of N standardized measurements such that 

a 0k ifijB lk^ a 2k a "-- Sz nk s0 

and 1„ 7 2 1 ^ „ 2 1 ^,2 1 ^ 2 1 

1T sZ 0k = Tr s 2 lk =_ sz 2k =...=~sz nk =l. 

We compute the correlation coefficients * 

r ij = r ji = lh Z ik Z jk (iandj = 0, 1 n) 

As a first approximation toward predicting z Q from z* , z„, . . . , z , we 
form the linear combination 



* Throughout the development the summation 2 will always be over k = 1 to N. 
From here on, however, the k subscript will be omitted for notational simplicity. 
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L = x lZl + x 2 z 2+ ... +x n z n , 

where x ., x,, . . . , x is any set of scalar constants. Later, the x become the varia- 
bles of simultaneous linear equations and should not be confused with the z which were 
first introduced as variables. The z are more precisely vectors; each z might, for 
example, be a set of test scores. In the final correlation analysis the z become im- 
plicit and the mathematics uses only their interrelations, the r... 

B. The correlation between z Q and L is derived as follows: 



i-2(z -m )(L-m L ) 



r 0L= 



where 



ff a L 



m o = I* sz o = ° 



a o =^f 2 ( z o- m o> 2 = 1 



so that r„ L jsj_ 



m L = TT SL = TT S(X 1 Z 1 + X 2 Z 2 + • •• + X n z n> 

= X l "i S Z l + *2 IT S Z 2 + ' • * + X n IT S Z r 
= Xl (0) + x 2 (0) + . . . + x n (0) 

= 

1 2z„L 



a L 



Tagging r QL with the new name MCC, dropping the subscript of ov, and 
letting 



= ¥ S Z ° L = ~N ^ Z (X 1 Z 1 + "2 Z 2 + " * + X n Z n) 



S 

N 



= x l -^ 2Z Z 1 + x 2^ Z Z 2 + ••• + *n ^ z z n 
= Xjr 10 + x 2 r 20 + . . . + x n r nQ 



provides the final relation 

MCC = 



S 
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C. It is now necessary to derive an expression for a or more conveniently, for 



a 2 : 



2 1 2 

a = 4r 2L (since m L = 0) 

= 1 , v2 

TT L ^ X 1 Z 1 + X 2 Z 2 + • • • + x n V 

= Xj ^ Zl (x lZl + x 2 z 2 + . . . + x n z n ) + . . . + x n isz n ( ¥l + . . . + x n z n ) 

= Xl (x^szj 2 + x 2 ^-S Zl z 2 + ... + x^SZjZ^h- ... + x n (x 1 ^ r Sz 1 z n+ ... + x n ^Sz n 2 ) 

= Xl (x 1+ x 2 r 12+ ... + Vln ) + -" +x n (x l r ln + *-' + V 



= x^a. + y^-n + + x_a. 



n n 



o 
The expanded representation of MCC or, more conveniently, of (MCC) is thus 

t \ 2 

(MCC) 2 = J 1 " 10 + *2 r 20 + ' ' • + x n r nO' 



x l a l + x 2 a 2 + '-- +X n a n 



D. An important question is which set or sets of x 1 , x„, . . . , x gives a maximum 
_ i & n 

(MCC) . This is answered by partial differentiation as follows: 



3 (MCC) 2 

9 Xj 


3 

ax. 


s 2 

a 2 - 


a 2 


a c 2 c 2 a 
a Xj s s S Xj 


■ a 2 




a 4 





But 3xy s2 = 5x7" (x l r 10 + • • ' + X n r n0) 2 

= 2Sr j0 
and 5 |_a 2 = 5 l-(x 1 a 1 + ... + x n a n ) 

= ( x ia!- a i + --- +X j3x-V"- +x n air a n > + a j 

Since &i = Xl r u + . . . + y.j + . . . + x n r. n , 

JL a -r 
9x7- a i r ij 

and 3^- ° 2 = < x l r jl + x 2 r j2 + " ' ' + Vjn> + a j 

= 2a. 

1 
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? 2Sa 2 r. n - 2S a. 
finally, g^-(MCC)^ = ] \ 1 

2 

If this is set equal to zero, it is seen that the minimum (MCC) occurs 
2 
when S, and hence (MCC) , equals zero, and the maximum when 

J-a-r.Q (j = l, 2, ..., n) 

But the least squares solution B- , B 2 , • . • , B is such that 

a i = r jo 

2 
This states that the solution obtained by maximizing (MCC) is just a constant mul- 

o 2 

tiple cr/S times the least squares solution. Now S and a are immediately available 

for any x 1 , x„, . . . , x . Therefore, once an x., x„, . . . , x is derived which maxi- 

2 
mizes (MCC) , the least squares solution can be quickly computed by the equation 

B i = X i ( 7 ) final 

2 

E. It remains only to outline the technique for converging upon the maximum (MCC) . 

One variable at a time is corrected upon. The correction is such that a maximal in- 

2 
crease in (MCC) is effected. This correction is derived as follows, assuming x p 

is the variable to be corrected upon: 

If 6 p is the correction such that the new x p = x p = x p + 6 p , 

then: the new S = S ; = x^a + • - • + ( x p + 6p) r pg + • • • + x n r n o 

= ( x i r io + • ■ • + x p r po + • • ■ + x n r n(P + 6 P r P0 
= S + 6pr p Q, 



in 



the new a. = a *. = x^^r .- + . . . + (x p + 6 p ) r. p + . . . + x n r. 

= < x l r il + ' ' • + x P r iP + • • * + Vin> + 6 P r iP 
= \ + 6 p r ip , 

and the new a =(cr ) = x.ai + . . . + (x p + 6 p ) a p + . . . + xa' 

= (x 1 a 1 + . . . + x R a n ) + 6 p (^r pl + ... + x^p^ + 6 p (a p + 6p) 



2 2 

= a + 6 p a p + 6 p a p + 6 p 

= or 2 + 2 ap 6 p+ 6 2 p 
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o/ 2 

Now, if (MCC) is the new (MCC) under a change 6 p in x p , 



(a 2 )^2s' J f-S')-(S') 2 3 §-(a 2 )' 

[ (* 2 )' f 

It should be emphasized that the partial differentiation is with respect to 6 p , the cor- 
rection in Xp to be determined. 
Using the fact that 

4;s'=4- (S + 6prpo, " rpo 

and 3 (o-V 3 ( <? + 2a p 6 p + 6 2 p ^ = 2 (a p + 6 p ) 

3 6p 35p 

we see that 

8 , Mpn 2' 2r po S^(cr 2 / -2(a p+ 6 p )(S') 2 
J3-CMCC) = 

[ (o 2 )' f 

Substituting in the derived expressions for S' and (a ) and setting this equal to zero 
gives the condition that 



g2r FQ - Sa p 
6 p = 



S - a p r po 

o 

Now 5 p is just the correction in x p that maximizes (MCC) with all other x. held 

constant. 

2 
As this process is repeated, each time correcting upon another x., (MCC) 

constantly increases and must approach the desired maximum. The method is there- 
fore convergent. 

The criterion for determining P at each stage is as follows: A 6 is calcula- 

2 
ted for each of the n x, using the S, cr , and the set of a that were computed after the 

prior correction. The largest numerical 6 is then selected as 5- and the corres- 

2 
ponding x p corrected upon by this amount. A new S and a is calculated, and on the 
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next iteration all e. are adjusted. 



e: = e. + 6„r. 



P'iP 



Since 



a i =r iO + e i ■ 



the new set of 6 can be computed simultaneously with the adjustment of the e.. 

2 2 1 

Besides an S and a a , an (MCC) is computed at the completion of each iter- 

o 

ation. Whenever the difference between successive (MCC) falls below an established 

tolerance, iterations are suspended. The x-, x„, .. ., x are obtained by summing 
the 6p, and they are multiplied by the ratio of the final S to the final <T to give 
Bj, B„, ..., B , the desired regression or least squares solution. 

It is interesting to note that when B*, B«, . . . , B n replaces the general 



n 



c l> x 2' 



x , the following relations exist: 



and 





2 

S= a = 


(MCC) 2 






a i =r iO 


(i= 1, 2, ..., n) 






e.= 






Table of Notation 




Possible 


Symbol 


Meaning 


Explanation 


Values 


pre 


proper r card 


which of three r cards to be 
read from 


0, 1, 2 


re 


r card 


number of r card going through 


0, 1, 2 


e# 


e code number 


setting up selectors on punch 




r# 


r code number 


setting up selectors on punch 




mr 


matrix row 


for e cards and r cards 


00 thru n-1 
for matrix of 
order N 


ef 


e field 


field of card from which e 
is read 


00 thru 24 


mc 


matrix column 


column to be read from on a 
given run 


00 thru n-1 


X ? e 


e card, with X„ 


exhausted e card; new e punched 

onXj 




X,Je 


e card, following 
X 7 


blank e card 





X 6 e 



e card with X 



6 



dummy e card 



X. designates that an 11 -punch appears in card column i. 
Y. designates that a 12-punch appears in card column i. 
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IV. Setup of the Decks 
F 5 Card 

The F6 card is the first card through the machine on each iteration and 
contains all of the data that is constant for calculations on the remaining cards. The 
F6 card contains the following data: 



Card Column 


Data 


1-2 


Variable identification 00 through n-1 


3-4 


e code number 


4 


12-punch to designate an F6 card 


5-9 


6 P 


10 - 14 


S 


15 - 19 


<r 2 


20 - 24 


(MCC) 2 


79 - 80 


iteration number 



The F6 card must be key-punched to start the problem. All succeeding 
F6 cards are punched by the 607. The variable having the largest r. n is corrected 
first, and consequently the initial F6 card appears as follows: 



Card Column 


Quantity 


Which is denoted 


1-2 


i 


mc 


3-4 


00 


ef 


4 


12 


identifies F6 


5-9 


r i0 


6 


10 - 14 


4 


S 


15 - 19 


4 


a 2 


20 - 24 


4 


(MCC) 2 


79 - 80 


00 


Iteration number 



It is important to remember that the initial i corresponds to the row which 
has the largest numerical r. Q . 

r Cards 



The r cards contain the elements r.. , with each card holding 25 three- 
digit elements of a row of the matrix. The present machine procedure provides for 
three r-cards or 75 variables. The r card contains the following data: 
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Card Column 


1 - 


2 


2 




3 - 


77 


78 - 


80 



Data 
Matrix row number, 00 through n-1 
12 -punch to indicate an r card 
25 r of a given row, 3 digits each 

r i0 



If n >25, more than one r card must be used for each row of the matrix. 
This requires the punching of more than one r deck. Twenty-five coefficients should 
be punched on the first deck, twenty-five on the second, and so on, until all have been 
punched. Since three positions are allowed for each r value, it is necessary to use 
.999 as the diagonal elements rather than 1.000. This has no appreciable effect on the 
final results. 

When more than one r card is required for each row, care should be exer- 
cised in adhering to the following codes for the various r cards: 



One r-card problem 
Two r-card problem 

Three r-card problem 



Second from 
Last r card 



X r 



Next to 
Last r card 



X, 



Last r card 
no additional punch 
no additional punch 

no additional punch 



e Cards 

The e cards are used for punching the successive errors on each iteration. 
A value e. is read from the ith e card and e^ is punched in the adjacent card field of the 
same card. All e cards have the following card form: 

Card Column Data 



1 
3 



2 

5 



matrix row number, 00 through n-1 
-r. (first e deck only) 

12-punch to identify as an e card 



There must be a "dummy" e card behind the F6 card. This "dummy" e 
contains an 11-punch in column 6 in addition to the 12-punch in column 3. This card is 
required to set up selectors for reading the first r card. 

The first set of e cards will contain the errors which result from first gues- 
ses of zeros for all weights. This results in errors e. = -r. , and these are punched 
in columns 3 - 5 as noted above. From then on the 607 reads e. from columns 3-5 and 
punches e( in columns 6-8, reads from 6-8 and punches in 9 - 11, and so on, until 
25 errors are punched on a card. On the 25th iteration (iterations are numbered 00-24) 

the machine punches an X 7 in the e cards and the machine stops. This indicates to the 

' / 

operator that the e cards are exhausted and that a new e deck of X „ must be inserted. 

This X „ deck is sorted or merged behind the X„ e deck. On the next run the e^ are 
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read from the last field of the exhausted e deck and the new e. are punched into col- 
umns 3 - 5 cf the new deck. The old e cards are then removed and the iterations are 
continued. When a third e deck is required, the exchange process is repeated. 

L6 Cards 

The L6 cards are blank cards containing the iteration number in columns 
79-80. One L6 card is punched at the end of each iteration, and an L6 card for the kth 
iteration becomes the F6 card for the (k + l)st iteration. 

B Cards 

These cards are blank and contain the variable number in columns 1-2. 
After the iterations have been completed, the final B.. , B«, . . . , B are punched in the 
B cards and are ready for listing. 

V. Running the Problem 

The final stack of cards, ready to be run through the machine, is in the 
following order face up: 

1. The first F5 card (number 00). 

2. A dummy e card containing an X fi in addition to the Y~. 

3. The r cards from the first row (row 00) of the matrix. Care must be taken to 
have the proper code punches as indicated previously. For a three r-card 
problem the first r card would contain an E ? , the second r card would contain 
an Xg, and the last r card would contain no code punches except the Y„ which 
identifies all r cards. 

4. The e card for the first row. 

5. Three and four above are repeated for all other rows. 

6. The L6 card. It will become the new F6 card. 

After every iteration a card moves from the stack of blank L6 cards to the 
bottom of the operational deck; the former L6 card becomes F6 and the former F6 card 
goes to the stack of used F6 cards. No control panel wiring changes are necessary 
throughout the entire operation. The machine automatically selects the appropriate 
field from each card, as directed by the r# and e#. 

When the e cards become exhausted, they will be replaced as previously 

described; the iterative procedure then continues normally. 

2 
The method will cause (MCC) to increase at each iteration, and it can be 

9 9 / 

stopped whenever the desired tolerance is achieved. As soon as I (MCC) - (MCC) I 
is within this tolerance, the machine will pick up calculator selector 13 and punch an 

X into column 1 of the L6 card. The operator may check the L6 and FS cards to de- 

2 2' 

termine if I (MCC) -(MCC) I is within the predetermined tolerance. 
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o 
The machine now places S/ct into the MQ unit. This constant is used to 

obtain B« , B , .... B„ from x 1 , x«, . . ., x It is necessary to inspect the MQ unit 

to see that it does contain the constant S/ct . If the MQ unit is blank, the quotient S/cr 

>1 , and the calculate panel must be changed to allow for one whole number and four 

decimals. This is done by making the following changes: 

Program 

102 Change Read Units Into 5 to Read Units Into 4. 

103 Change Read Units Out of 6 to Read Units Out of 5. 
137 Change Read Units Into 6 to Read Units Into 5. 

After the three wire changes are made, the punch panel is opened to drop 

out calculator selector 13 and the last iteration is repeated to recompute the constant 

9 2 

S/ct . As soon as calculator selector 13 is picked up and S/ct is located in the MQ 

unit, the least squares weights may be computed. 

In order to obtain the final desired weights, it is necessary to sort the 6 
cards in variable number order with a blank B card behind each of the corresponding 
group of 6 cards. This deck is then run through the 607, which has been set up by 
calculator selector 13, to compute the sum of 6 corrections for each variable and mul- 
tiply by S/ct 2 to obtain the B weights. The B weights are then removed from the deck 
and are ready for listing. It may be desired to check computations by substituting the 
B weights into the equations. This can be done by using the B weights as a predeter- 
mined set of F6 cards and placing a blank r card (Y 2 card) at the end to keep calcula- 
tor selector 13 from transferring. Thus the checking will require as many iterations 
as there are non-zero B weights. At the end of the last iteration, the last field that is 
punched in each e card should be approximately equal to zero. 

VI. Description of Programming 

Section VI comprises a detailed program of the procedure for the Type 607, 
together with auxiliary charts and tables. It is too bulky for inclusion in this article, 
but will be made available upon request. * 

VII. Solution of a 24th Order Regression Weight Problem 

As was pointed out earlier, up to 75 variables can be treated by the program 
without further modification. With some changes this upper limit may be raised arbit- 
rarily high. In general, convergence will be slowed down somewhat as the size of the 



* Requests for the programming description should be directed to: 

The Applied Science Division, 
International Business Machines Corporation, 
590 Madison Avenue, New York 22, N. Y. 
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problem increases. 

To provide an indication as to the convergence time of this method, a rep- 
resentative problem of the 24th order has been selected for discussion. Thirty-two 
iterations were required for convergence to an MCC, constant to four places. When the 
derived set of B values were substituted back in the given equations, the magnitude of 
the largest error was .007. The overall running time was under an hour. In general, 
the amount of time required is a direct function of the number of iterations necessary 
for the desired accuracy. Since the MCC is computed as part of the machine procedure 
after each iteration, an immediate measure of convergence is always automatically 
available. 
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AUTOMATIC DIAGNOSIS OF CPC CODING ERRORS 

D. W. Marquardt 

E. I. du Pont de Nemours and Company, Inc. 
Wilmington, Delaware 

Introduction 

One of the major costs in the use of digital computers for technical computa- 
tion is that of "debugging" new programs. Accordingly, the development of systema- 
tic procedures for the diagnosis of coding errors offers the possibility of significant- 
ly reducing computation costs. A description will be given here of some techniques 
that have been in routine use with marked success for more than a year. These tech- 
niques, although designed for use with an adaptation of the Heising floating decimal 
setup, are applicable with only minor modifications to many of the general purpose 
setups currently used with the Model II CPC. These setups typically provide for up to 
three operands, leading to a single result, on each card cycle. 

Clas si fication of Errors 

In general, there are two broad classes of errors that may arise in the prep- 
aration of problems for a digital computer. First, there are logical errors in plan- 
ning the strategy of problem solution. Errors of logic are not discussed here. 

Secondly, there are coding errors which arise in the transcription of the de- 
tailed strategy of solution into computer code. CPC coding errors may be grouped con- 
veniently into two types as follows: 

1. Type 1 Errors 

These errors result from violations of built-in limitations of the CPC 

(with a typical setup), namely: 

(a) The addresses specified for the three operands on any card must 
reflect the limitations that not more than two operands can come 
from 941 -type storage, ncr mere than two from accounting mach- 
ine storage. 

(b) The addresses specified for the three operands and for the result 
must satisfy the timing restrictions, which are as follows: 

(1) It is not possible to read an operand out of an accounting 
machine storage location and read the result back into that 
same location. It is possible, however, to do this with 
941 -type storage locations. 
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(2) The result of the (n - l)st card, if sent to an accounting 
machine storage location, may not be obtained from that 
location as an operand on the nth card. 

(3) The result of either the (n - 2)nd or (n - l)st card, if sent 
to a 941 -type storage location, may not be obtained from 
that location as an operand on the nth card. 

(4) The results from two successive cards may not be sent 
to the same storage location. 

2. Type 2 Errors 

This category includes all coding errors not due to the limitations listed 
above . 

Significance of Type 1 Errors 

A major significance of Type 1 errors is that they are readily susceptible to 
automatic diagnosis. Their automatic diagnosis is made more important by the fact 
that some Type 1 errors will not always, or even most of the time, cause the CPC to 
malfunction. Consequently, they will frequently escape detection during normal de- 
bugging of a code but manifest themselves during the computational phase. The re- 
mainder of this article will describe techniques for the automatic diagnosis of Type 1 
errors. 

Techniques of Diagnosis 

It is possible to check for the presence of coding errors by means of a spec- 
ial run or runs (on the CPC or on some other machine*) prior to the computation runs. 
It is alternatively possible to have the CPC perform at least some of these checks con- 
tinuously as the computations proceed. 

There is much to recommend the latter procedure. The CPC is frequently 
employed for engineering-type calculations. These are often characterized by the self- 
service use of the CPC by non- specialists, and by changes in problem complexion as 
the results develop, necessitating changes in strategy and hence changes in coding. Ex- 
perience has shown that coding errors are somewhat more frequent in coding modifica- 
tions than in original coding. Therefore, a continuous check is desirable, provided it 
is not necessary to sacrifice a significant part of computing capacity in order to obtain 
it. 

The present procedure is a compromise between prechecking and continuous 
checking. A continuous check is made for Type 1-a errors at every card cycle. 

*Note: The present procedures employ the CPC exclusively. A control panel to per- 
form some of the checks described here has been developed using the Type 101 Statis- 
tical Machine (1). Seme of the checks described here have also been included in the 
Rand Corporation "Duplex" setup (2). 
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A special 605 panel and a 527 panel to go with it are used to make a precheck 
for all timing errors (Type 1-b) in a single pass of the instruction cards through the 
527 punch unit. 

It will be assumed in what follows that not more than 6 banks of 941 storage 
are present and that accounting machine storage has addresses 91, 92, . . . , 98. The 
manner in which the system should be modified for variations of this will be obvious. 

Control Panel Wiring to Check for Type 1-a Errors 

The accounting machine panel wiring to perform this check is shown in Figure 
1. The CPC will stop and the stop light will turn on if any card fails to meet the re- 
quirement that there be at least one "9" but not more than two "9's" in the tens posi- 
tions of the addresses for the three operands. Obviously, it will be necessary to punch 
a dummy "9" in any blank cards. To restart the CPC, the final total button and then 
the start button must be pressed. 

It will be noted that relatively few components are used to perform the check, 
so that there is no significant loss in computing capacity. 

Control Panels to Check for Type 1-b Errors 

These panels will automatically check a deck of instruction cards for the pre- 
sence of all Type 1-b errors (timing restriction violations). It is recommended that 
the cards be passed through this timing check procedure before debugging is begun, 
and again after it is otherwise completed. The panels (527 and 605) require only the 
605 unit and the 527 unit of the CPC. The cards are passed through the 527 unit, and 
if no timing errors are present, the cards will not stop. However, if a timing error 
is present, the zero check light on the 527 unit will light, and subsequently the cards 
will cease to feed. The last card to be ejected into the stacker before the cards stop 
feeding contains the timing error.* If an error has been detected, so that the zero 
check light is lit, then the reset button must be pressed before the start button can 
cause feeding to resume. 

To clear the 605, the start button must be depressed for several cycles before 
the cards are placed in the feed hopper. It may happen that the zero check light (some- 
times also the unfinished program light) will turn on during these clearing cycles. If 
so, then the reset button must be depressed following each cycle. 

The 605 panel wiring is given in Figure 2, and the 527 panel wiring in Figure 
3. Figure 4 shows the flow chart for the 605 panel wiring. 

♦Note: If an error is detected near the end of a deck, so that the card-feed hopper emp- 
ties and causes feeding to stop before the error card is ejected into the stacker, then 
the start button should be pressed until the zero check light causes feeding to stop, with 
the error card the last one ejected in the stacker. 

This eventuality may be avoided, if desired, merely by placing two blank cards 
at the end of the deck of instruction cards. 
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Figure 2 



STEP 

1 
2 
3 
4 

5 

6 

7 

8 

9 
10 
11 
12 
13 
11 

15 
16 
17 
18 

19 
20 
21 
22 
23 
24 
25 
26 
27 
28 
29 
30 
31 
32 
33 
34 
35 
36 
37 
38 
39 
40 
41 
42 
43 
44 
45 
46 
47 
48 
49 
50 

51 
52 
53 
51 
55 
56 
57 
58 



SUPPRESSION 



N.Z. 



N.Z. 



N.Z. 



N.Z. 



N.Z. 



N.Z. 



N.Z. 



N.Z 



N.Z. 



N.Z. 



N.Z. 



N.Z. 



(-) 



(1) (2) 





FUNCTIONS 




1. 


2. 


3. 


Emit 8 


CRI - 


Into 2 


F4R0 


CRI + 


B.T. 


Reset 






F1R0 


CRI + 




F4R0 


CRI - 




Z.T. 


Reset 






(P.U. GR. 


SUPP l) 


P2R0 


CRI + 




F4R0 


CRI - 




Z.T. 


Reset 






(P.U. GR. 


SUPP l) 


F3R0 
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F4R0 


CRI - 




Z..T. 


Reset 
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Z.T. 


Reset 
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CRI - 




Z.T. 


Reset 






(P.U. GR. 


SUPP 2) 


F3R0 


CRI + 




G2R0 


CRI - 




Z.T. 


Reset 






(P.U. GR. 


SUPP 2) 


G2R0 


CRI + 




Z.T. 


Reset 






(D.O. GR. 


SUPP 2) 


GIRO 


CRI + 




Z.T. 


Reset 






(P.U. GR. 


SUPP 4) 


F1R0 


CRI + 




GIRO 


CRI - 




Z.T. 


Reset 






(P.U. GR. 


SUPP 1) 


F2R0 


CRI + 




GIRO 


CRI - 




Z.T. 


Reset 


SUPP 1) 




(P.U. GR. 


F3R0 


CRI + 




GIRO 


CRI - 




Z.T. 


Reset 






(P.U. GR. 


SUPP 1) 


GIRO 


CRI + 




F4R0 


CRI - 




Z.T. 


Reset 






(P.U. GR. 


SUPP 1) 


G2R0 


CRI + 




Emit 8 


CRI - 


Into 2 


B.T. 








(D.O. GR. 


SUPP 2) 


GIRO 


G2RI 




F4R0 


G1RI 


Reset 


Emit 3 


CRI - 




Emit 5 


CRI + 




(B.T. FOR 


SEL. P.U.) 




605 PANEL WIRING FOR 






TIMING CHECK PROCEDURE 
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Figure 4a 
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FLOWCHART -TIMING CHECK PANELS 
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ERRATA 
Technical Newsletter No. 8 



Page 13: Paragraph 4 should read: At a card punching rate of 100 cards 

per minute, up to 554 ms are available for computing; and at a card reading speed 
of 200 cards per minute, up to 254 ms are available for computing. (The interlock 
timing for the 533 is 46 ms for both read and punch. ) 

Page 22: It has been discovered that the following locations are not used 
by the floating decimal interpretative system: 0047, 0173, 0319, 0353, 0372, and 
0423. Locations 0019-0024 are not used either, but they were purposely left blank 
in order to allow the addition of 6 more operations if so desired. Location 0466 is 
used by the system, however, making a total of 467 memory locations required. 

Page 23: General interpretation, line 12 should read: 

0040 SU 11 0193 8003 

Page 24: Problem description at the top of the page should read: 

(A) + (K) > K Subroutine 

(K) = a. A, , (A) = a„ A„, a = xx, A = x.xxxxxxx± 

Page 25: In the floating decimal multiplication subroutine, instruction 0181 
will result in an overflow in the case of all 9*s. This will in turn invalidate the two 
subsequent Store Address instructions. 

Page 29: Interpretation of 07, lines 14 and 15 should read: 

0243 RAL 65 B 0465 

0465 STL 20 0089 0134 

Page 31: Interpretation of 11, line 2 should be replaced by the following 
2 instructions: 

0417 RAL 65 8002 0425 

0425 AU 10 0037 0378 

Page 31: Interpretation of 13,lines 1 and 2 should be replaced by the 
following 3 instructions: 



0013 


SRT 


30 


0002 


0271 


0271 


RAL 


65 


8002 


0229 


0229 


AU 


10 


0037 


0431 



Page 34: Interpretation of 16, line 35 should read: 

0326 2100000026 

Page 36: Interpretation of 18, lines 1-4 should be replaced by the following 
3 instructions: 

0018 SRT 30 0002 0384 

0384 STL 20 0389 0442 

0442 RAL 65 0029 0401 

Page 40: First line should read: Store next instruction in 0083. 

Page 41: In the Remarks column, a, = a, = a„ -1 should be replaced by 
a 3 = a 2 " a l " L 

Page 61: It has been discovered that location 0065 is not used by the double 
precision arithmetic program. 
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Technical Newsl et ter No. 9 
Page 27: 



ERRATA 



Line 9 should read: A control word of the form n n + A Q A Q . 



Page 34: 
sine subroutine. 



1.6., = for the sine subroutine should read: 1.0J = 6 for the 



Page 35: It has been discovered that location 0077 is not used by the sine- 
cosine subroutine SCI. 

Page 37- Line 3, opposite instruction location 0006 should read: Compute 
6 = tt/2 - 8. If 6 n 10, subtract 2ir from it. Note: For the cosine subroutines, it 
is necessary to place -6 in the upper half of the accumulator. 

Page 51: Paragraph 2, lines 1-2 should read: This subroutine is initiated 
by placing the argument x= x. xxxxxxxxx in the upper accumulator. 

It has been discovered that location 0023 is not used by the logarithm sub- 
routine. 

The Matrix- Vector multiplication sub-routine requires that all factors be less 
than 1 in absolute value. Thus a.., x. = .xxxxxxxxxx 

The following table lists the required contents of the accumulator when entering 
the various sub-routines and the results left in it upon their completion. 



Sub- Routine 

SRI Ei 

E: 


Upper 
try 


Lower 
A 


it 


VHT 


SR2 


Can be anything 


A 





VTT 


SCI, SC2 


e 








sine 8 or cos 8 


Arc sin 


X 


Can be anything 


arc sin x 


intermediate results 


Arctan 


X 





arctan x 


intermediate results 


Arctan by 
Radix Table 





X 


arctan x 


intermediate results 


e x 


X 


Can be anything 





e* 


In (1+x) 


X 





In (1+x) 


intermediate results 
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